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1.  INTRODUCTION 


The  Air  Force  Geophysics  Laboratory  (AFGL)  has  been  tasked  by  the 
Space  Division  (SD)  of  the  Air  Force  Systems  Command  to  provide  meteoro¬ 
logical  support  and  analyses  in  support  of  the  Advanced  Ballistic  Reentry 
Systems  (ABRES)  Program.  The  primary  objectives  of  the  work  carried  out 
under  this  contract  by  Environmental  Research  $  Technology,  Inc.  (ERT) 
have  been  to  assist  AFGL  in  (1)  the  design  and  analyses  of  experiments 
to  assess  the  probable  error  of  the  AFGL  aircraft/radar  methodology  for 
characterizing  the  cloud  microphysical  environment  along  a  reentry 
trajectory,  (21  the  development  and  testing  of  improved  analytical 
methods  for  des<  ribing  the  cloud  microphysical  environment  along  a 
reentry  trajectory,  and  (3)  the  generation  of  environmental  statistics  to 
be  used  as  input  to  operations  analyses  of  possible  future  reentry  tests 
of  experiments.  Another  major  portion  of  the  work  involved  the  mainten¬ 
ance  and  running  of  many  computer  programs;  these  programs  have  been 
documented  an  I  form  a  part  of  this  report. 

Throughout  the  two-year  period  of  the  contract,  the  work  was  broken 
down  into  six  tasks.  The  task  titles  are: 

1)  Cloud  and  Precipitation  Descriptions  for  Ten  Case  Studies, 

2)  Assessment  of  Radar  Parameters  for  Hydrometeor  Scatter 
Measurement  Using  the  ALCOR  and  TRADEX  Radars  at 
Kwajalein, 

3)  The  Weather  Correlation  Runs  of  June  26,  1979  at 
Kwajalein, 

4)  Cirrus  Cloud  Particle  Detectability  Using  the  Kwajalein 
Radars , 

5)  Assessment  of  Airborne  Radar  AN/APQ-1 22 (V)  8  for  Use  in 
Quantitative  Measurements  of  Liquid  Water  Content,  and 

6)  The  Processing  of  Radar  and  Aircraft  Weather  Data  from 
Wallops  Island. 

This  report  describes  the  key  results  of  each  of  the  six  tasks.  In  addi¬ 
tion,  the  procedure  for  assigning  the  liquid  water  content  in  various 
atmospheric  layers  is  presented  in  Appendix  D. 
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2.  ANALYSES  AND  SUPPORT 


2.1  Task  1  -  The  Processing  of  Radar  and  Aircraft  Weather 
Data  from  Wallops  Island 

An  important  part  of  this  program  was  the  maintenance  and  running 
of  many  different  types  of  computer  programs.  The  programs  were 
designed  to  process  weather  radar  and  aircraft  data  of  clouds  and  pre¬ 
cipitation.  There  were  10  major  programs  which  were  run  at  various 
intervals  over  the  two-year  span  of  the  program,  and  a  total  of  more 
than  500  computer  runs  were  made. 

Documentation  of  the  computer  programs  used  throughout  the  program 
is  given  in  Appendix  A.  A  listing  of  all  of  the  computer  programs  is 
available  from  the  Meteorology  Division  of  the  Air  Force  Geophysics 
Laboratory. 

2.2  Task  2  -  Assessment  of  Radar  Parameters  for  Hydrometeor 
Scatter  Measurements  Using  the  ALCOR  and  TRADEX  Radars  at 
Kwajalein 

The  results  of  this  task  have  been  described  in  a  report  by  Blood 
and  Crane  (1979).  The  report  documents  the  scattering  relationships  and 
radar  parameters  which  enter  into  the  measurement  of  effective  radar 
reflectivity,  Ze>  obtained  from  the  ALCOR  (C-band)  and  TRADEX  (S-band) 
radars  operated  at  Kwajalein  Atoll  for  hydrometeor  scatter  investigations. 
An  assessment  is  given  of  the  calibration  procedures  used  for  the  radars 
including  assumptions,  corrections,  and  uncertainties  associated  with 
the  parameters  used  in  the  measurements.  The  current  best  estimate  radar 
constants  for  ice  and  water  used  for  data  processing  are  documented  for 
users  of  the  data  together  with  an  editing  criteria  for  rain  data  based 
upon  the  principal -to-opposite  polarization  received  power  ratios. 

The  summary  and  conclusions  of  the  report  are  as  follows: 

"The  ALCOR  (C-band)  and  TRADEX  (S-band)  radars  at  Kwaja- 
lcin  are  well-calibrated  instrumentation  radars  which  serve 
among  many  other  roles  to  provide  measurements  of  the  environ¬ 
mental  conditions  of  hydrometeor  activity  within  the  reentry 
corridor  of  the  Kwajalein  Missile  Range  (KMR) .  The  report 
defines  the  essential  radar  parameters,  assesses  significant 
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sources  of  errors  in  these  parameters  and  documents  the  most 
recent  best  estimate  of  the  appropriate  radar  constants  to 
use  in  data  processing.  A  discussion  of  the  calibration  pro¬ 
cedures  is  also  given  with  a  recognition  of  the  ultimate 
importance  of  sphere  calibrations  in  close  space  and  time 
proximity  to  the  hydrometeor  scatter  measurements  of  interest 
for  minimizing  the  calibration  uncertainties.  The  overall 
RMS  (lo)  uncertainties  associated  with  the  radar  measurement 
are  estimated  at  1.0  dB  each  for  ALCOR  and  TRADKX.  This 
assumes  that  the  bias  terms  associated  with  signal  process¬ 
ing  are  appropriately  accounted  for  by  the  use  of  the  best 
estimate  correction  factor  and  that  the  radar  parameters  do 
not  change.  Implicit  in  this  analysis  is  the  assumption  that 
a  zero  dB  polarization  mis-match  loss  occurs  for  hydrometeor 
scatter  and  that  excess  losses,  such  as  radome  attenuation 
changes  due  to  rain  or  condensation  of  moisture  on  the  ALCOR 
radome,  do  not  occur  during  the  measurements. 

An  initial  attempt  to  classify  rain  scatter  on  the 
basis  of  principal-to-opposite  polarization  ratio  is  presented. 

This  tool,  useful  for  data  editing  and  interpretation,  sug¬ 
gests  that  further  means  of  classifying  the  hydrometeor  back- 
scatter  may  be  intrinsically  available  in  the  different  radar 
polarizations  and  worthwhile  for  further  investigation  at 
Kwajalein." 

2.3  Task  3  -  The  Weather  Correlation  Runs  of  June  26,  1979 
at  Kwajalein 

The  purpose  of  the  1 inked-radar/aircraft  experiment  on  June  26, 

1979  at  Kwajalein  was  to  collect  data  simultaneously  by  radar  backscatter 
and  by  on-board  aircraft  measurements  of  the  hydrometeor  type  and  size 
distribution  for  correlation  and  calibration  purposes.  The  details, 
results,  and  recommendations  of  the  experiment  are  provided  in  Appendix  B 

2.4  Task  4  -  Cirrus  Cloud  Particle  Detectability  Using  the 
Kwajalein  Radars 

The  ability  to  make  measurements  of  low  particle  densities  of 
tenuous  clouds  by  radar  is  of  interest  to  the  question  of  reentry 
vehicle  surface  erosion.  A  report  on  the  detection  of  low  concentration 
particles  within  cirrus  clouds  is  given  in  Appendix  C.  The  conclusion 
of  the  study  is  that,  using  the  multi- frequency  radars  at  Kwajalein, 

-3 

80  pm  ice  particles  with  densities  as  low  as  42  m  could  be  detected  at 
a  range  of  42  km  using  the  TRADF.X  L-band  chirp  system. 
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2.5  Task  5  -  Assessment  of  Airborne  Radar  AN/APQ-122  V)  8 

for  Use  in  Quantitative  Measurements  of  Liquid  Water  Content 

2.5.1  Introduction 

This  report  summarizes  the  findings  of  an  investigation  into  the 
potential  scientific  use  of  a  multi-function  airborne  radar  to  be 
referred  to  as  the  APQ-122,  which  is  currently  operating  on-board  an 
instrumented  C-130  aircraft.  The  radar  is  being  considered  for  use  in 
quantitative  measurements  of  liquid  water  content  in  addition  to  its 
currently  designed  modes  of  weather  mapping  and  other  functions.  In 
order  to  perform  these  data  collection  measurements,  an  investigation 
was  initiated  toward  possibly  using  the  radar  outputs  for  making  reflec¬ 
tivity  measurements  in  a  pulse  volume  just  ahead  of  the  aircraft.  The 
purpose  was  to  define  the  probable  use  of  the  instrument  without  major 
modification  of  existing  hardware  and  to  define  auxiliary  sampling  and 
recording  equipment  that  would  be  needed  for  collecting  the  desired 
measurements. 

An  initial  look  at  the  sensitivity  of  the  radar  at  Ka-  and  X-bands 
(Blood,  1980j  seemed  to  indicate  that  the  desired  levels  of  liquid  water 
content  would  be  detectable  providing  that  the  radar  output  could  be 
sampled  at  ranges  of  nominally  1-2  Kft  ahead  of  the  aircraft.  Mr.  Vernon 
(5.  Plank  (private  communication)  indicated  that  detectable  thresholds  of 
0.003  g/m'^  at  Ka-band  and  0.001  g/m'>  at  X-band  (long  pulse)  should  be 
adequate  for  measuring  cirrus  type  particles. 

Further  information,  however,  was  needed  to  assess  the  performance 
to  be  expected  from  the  instrument.  Due  to  the  lack  of  immediate  access 
to  equipment  manuals,  specifications  and  any  available  specific  measure¬ 
ments  of  the  type  needed  to  completely  describe  the  APQ-122  radar  per¬ 
formance  at  near  range,  ERT  requested  that  additional  information  be 
obtained  through  the  Radar  Branch  at  Wright  Patterson  AFB.  The  AF 
responses  and  flight  manual  specifications  have  been  used  to  help  specify 
additional  tests  that  may  be  needed  and  to  propose  a  tentative  plan  of 
radar  instrument  usage,  should  all  of  the  assumptions  be  verified.  It 
currently  looks  promising  for  the  quantitative  data  to  be  collected; 
however,  several  unknowns  about  performance  remain.  The  questions 
that  need  further  verification  are:  (1)  are  the  Ka-band  and  X-band 
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transmitter/receiver  systems  stable?  (2)  does  the  radar  have  adequate 
sensitivity  and  dynamic  range  at  the  near  ranges  (clutter  free)?  and  (3) 
can  the  signal  amplitude  be  calibrated  with  minimal  bias  errors  and 
uncertainties  (M  to  2  dB  each)? 

With  the  APQ-122  operating  in  the  weather  mode  (pencil  beam)  and  the 
antenna  in  a  stop  mode  oriented  ahead  of  the  aircraft  and  compensated  for 
aircraft  drift  angle  on  a  horizontal  flight  path,  current  findings  are 
that  the  radar  should  be  able  to  sense  the  reflectivity  levels  (ampli¬ 
tude  only)  at  ranges  between  270  m  (900  ft)  and  750  m  (2500  ft).  The 
APQ-122  contains  two  completely  separate  transmitter/receiver  units.  As 
described  in  the  operator's  flight  manual  LTM1MC-1301;-1 ,  both  the  Ka-band 
(v33  GHZ)  and  X-band  (a.9.36  GHz)  radar  frequencies  may  not  be  operated 
simultaneously  as  the  mapping  antenna  is  shared  in  the  weather  mode; 
thus  the  X-band  mapping  functions  are  not  available  during  Ka-band 
mapping;  i.e.,  the  Ka-band  mapping  subsystem,  when  operating,  has 
priority  over  the  X-band  system,  and  X-band  automatically  reverts  to 
standby.  The  pencil  beam  mode  is  selectable  (versus  fan  beam),  and  the 
antenna  servos  receive  positioning  information  from  the  electronic 
control  amplifier  in  the  form  of  tilt,  roll,  pitch,  drift,  and  azimuth 
drive  commands  which  compensate  for  airplane  maneuvers.  The  tilt  (verti¬ 
cal)  angle  is  selectable  from  10°  upward  to  15°  downward  from  horizontal 
with  the  aircraft. 

The  exact  use  of  the  two  (Ka-  and  X-)  radar  bands  has  not  been 
completely  determined,  and  usage  may  depend  on  the  type  of  particles 
being  sampled.  It  may  be  desirable  to  use  either  radar  depending  on  the 
altitude  being  studied.  The  33  GHz  frequency  may  be  better  for  ice 
measurements  at  higher  altitude  and  the  9  GHz  frequency  more  optimum  for 
rain  at  lower  altitudes.  At  Ka-band  the  reflectivity  of  ice  is  higher 
than  at  X-band;  however,  when  within  rain,  Ka  frequency  is  not  as 
desirable  as  X-band  due  to  increased  wave  propagation  attenuation.  The 
option  of  using  either  system  may  be  desirable  pending  experimental 
observations  at  which  time  an  operational  procedure  may  be  established. 

The  following  sections  describe  the  radar  sensitivities  expected, 
preparatory  test  recommendations  (validation) ,  and  a  tentative  outline 
of  data  collection  requirements.  Attempts  have  been  made  to  indicate 
the  assumptions  where  quantitative  measurements  are  not  yet  available. 


In  order  to  guarantee  a  successful  experiment,  pric r  ro  the  fabrication 
of  the  hardware  each  uncertainty  should  be  checked  either  by  bench 
measurement,  within  the  aircraft,  or  against  such  manufacturers  tests 
which  may  already  exist  on  the  specific  unit  to  be  flown. 

2.5.2  System  Sensitivity  to  Detect  Ice  and  Water 

The  Ka-band  radar  (0.2  ns  pulse  width)  sensitivity  was  estimated  by 
Blood  (1980);  it  was  based  on  quoted  nominal  values  of  transmitted  power, 
antenna  beamwidths  and  receiver  sensitivity  at  pulse  repetition  frequen¬ 
cies  (PRFs)  of  2000  and  4000  PPS.  An  integration  gain  of  17  dB  (0.589 
sec)  was  obtained  at  4000  PPS  whereby  2560  pulses  would  be  averaged 
incoherently  during  the  time  that  2.65  resolution  volumes  would  be  swept 
out  by  the  aircraft  at  a  260  Knot  indicated  air  speed  (IAS).  The  result¬ 
ing  level  of  dRZe  detectable  for  ice  at  a  range  of  1000  meters  was 
-IS  dBZe  at  a  10  dB  S/N  level.  For  water,  a  6  dB  lower  dBZe  level  was 
obtained  under  the  same  conditions.  If  the  lower  (2000  PPS)  PRF  were  to 
be  used,  the  processing  gain  proportional  to  Jx  would  drop  to  15.5  dB, 
and  detectable  levels  would  be  reduced  by  1.5  dB.  In  considering  the 
system  losses  to  arrive  at  these  values,  a  nominal  5  dB  value  was  esti¬ 
mated  (Blood,  1980).  This  loss,  typical  for  airborne  weather  radar, 
usually  breaks  down  as  follows:  radome,  1.5  dB;  waveguide  and  mismatch 
losses,  1.0  dB;  antenna  losses,  1  dB;  and  beam  shape  loss,  1.5  dB.  The 
actual  system  losses  of  the  APQ-122,  however,  are  currently  unknown. 
Quoted  values  were  used  for  other  parameters  including  a  receiver  minimum 
discernable  signal  (MDS)  of  -90  dBM.  Not  considered  in  the  receiver 
sensitivity  however  was  the  existence  of  a  built-in  (fixed)  STC,  sensi- 
tivitv-timc  control  taper.  The  STC  depth  was  assumed  to  be  manually 
controllable  in  the  weather  mode  and  therefore  settable  to  tero.  Since 
the  STC  values  were  completely  unspecified,  nothing  can  be  definitely 
stated  about  the  K-band  radar's  capability  in  achieving  the  above  reflec¬ 
tivities  based  on  the  -90  dBM  MDS  value.  Further,  since  no  actual  system 
losses  were  given,  the  departures  from  the  estimated  5  dB  value  may  be 
several  dB. 

At  X-band,  the  detectability  was  based  on  the  most  sensitivity 
(long  pulse)  waveform  of  4  jjs  duration.  In  view  of  the  probable  receiver 
recovery  estimate  of  at  least  4  pulse  widths  (i.e.,  16  ,,s  =  2400  m) , 


6 


this  was  perhaps  an  unrealistic  waveform  choice  for  X-baiul.  The  narrower 
pulse  width  (0.6  and  0.3  ys)  X-hand  waveforms  have  therefore  been 
reexamined  in  Table  1  along  with  those  previously  reported  at  Ka-band. 
Choosing  the  higher  PRF  (narrowest  pulse)  waveform,  the  integration  gain 
for  4096  pulses  is  about  18  dB.  During  this  integration  time,  6.1 
resolution  volumes  are  swept  out.  With  an  18  dB  processing  gain,  a 
revised  X-band  sensitivity  estimate  is  given  in  Table  2  (lower  portion) 
for  the  0.3  s  pulse.  The  X-band  radar  is  less  sensitive  with  this 
waveform  than  at  4  us  (MDS  =  -98  versus  -106  dBM)  and  consequently, 
because  of  the  lower  sensitivity  and  scattering  volume  size,  the  level 
of  reflectivity  is  reduced  considerably  (-10  dBZe  versus  -25  dBZe  for  ice 
at  4  y s  and  1000  M) . 

In  all  probability,  the  receiver  outputs  at  Ka-  and  X-bands,  as 
mentioned  previously,  may  be  sampled  at  ranges  between  270  m  and  750  m. 
Using  the  750-meter  value  as  a  conservative  estimate,  the  detectability 
at  Ka-band  for  ice  would  be  -20.5  dBZe  ('v. 003  g/m‘>  LIVC)  and  at  X-band, 
-12.5  dBZe  (%.008  g/m'^  LWC).  A  more  exact  definition  of  the  radar 
detectability  may  be  made  once  the  near-range  receiver  response  charac¬ 
teristic  has  been  determined  from  measurement.  Tt  is  necessary  that  all 
feed-through  modulation  and  near-range  "clutter"  type  returns  have 
decayed  to  the  receiver  noise  level  quoted  or  the  sensitivities  will  be 
degraded  correspondingly.  The  radar  parameters  used  in  these  detecta¬ 
bility  estimates  are  given  in  Table  3  including  the  site  of  the  illumi¬ 
nated  volume  assuming  a  filled  beam  and  single  pulse  width  resolution  at 
a  range  of  "750  meters  ahead  of  the  aircraft.  Blood  (1980)  describes  all 
of  the  relationships  used. 

2.5.3  Preparatory  Test  Recommendations 

Transni tter/Recei ver  Stability 

I.iuu'. ■!  w.i*er  content  measurements  are  derived  from  effective  radar 
reflect iv:ty  measurements  (dBZe)  which  in  turn  are  derived  from  calibrated 
amplitude  -‘'asurerents  of  the  received  power  from  a  gated  portion  of  the 
’-eceiver  out  nut.  These  measurements  therefore  depend  on  amplitude 
stability  of  the  transmitter-receiver  system  both  pulse-to-pulse  and 
over  a  long  period  of  time.  It  is  probable  that  magnetron  power  outputs 
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will  vary  or  fluctuate  especially  if  impedance  variations  or  temperature 
changes  within  the  feed  and  antenna  system  occur.  The  characteristics 
of  the  transmitter  system  may  be  monitored  partially  by  bench  tests  ahead 
of  time  and  by  in-flight  system  measurements.  The  receiver  gain  stability 
needs  also  to  be  known.  This  should  be  obtainable  mainly  from  bench 
tests  if  the  data  do  not  already  exist. 

Sampling  Range  and  Recovery  Time 

It  is  desired  to  know  the  video  channel  time  response  in  the 
vicinity  of  the  transmitted  pulse  in  order  to  determine  the  optimum 
gating  point  for  L.W.C.  measurements.  Several  effects  may  determine 
this  gating  point.  The  video  is  blanked  to  1.8  ps  (in  the  case  of 
X-band)  beyond  the  start  of  the  pre-master  trigger  (video  blanking)  and 
modulator  trigger  pulses.  The  receiver  output  therefore  cannot  be  gated 
prior  to  270  meters  range.  Beyond  this  point,  built  in  STC  circuits 
affect  the  receiver  gain  and  ringing  responses  may  sti 1 1  occur  from  the 
trailing  edge  of  the  transmitted  pulse  or  modulator  feed  through  signal. 
Any  other  "clutter"  returns  occurring  in  the  operating  system  output 
during  this  time  are  usually  by  multiple  reflections  in  the  feed,  antenna 
system,  or  at  any  impedance  discontinuities  experienced  by  the  RF  energy. 
Therefore,  the  shape  of  the  video  output  must  be  measured  quantitatively. 
Measurements  could  be  measured  on  an  A-scope  display  (wide  bandwidth 
oscilloscope)  to  determine  the  response  shape  and  desired  sample  gate 
position  where  returns  are  sufficiently  close  to  the  noise. 

Receiver  Dynamic  Range  and  Sensitivity 


Once  the  sample  gate  positions  have  been  obtained  for  Ka-  and  X-bands, 
the  receiver  dynamic  range  and  sensitivity  may  be  measured.  The 
receivers  have  logarithmic  IF  responses  and  no  automatic  gain  control 
(i.e.,  fixed  gain)  circuitry.  Therefore,  with  a  pulsed  calibrated  RF 
signal  generator  together  with  a  wide-bandwidth-DC-coupled  calibrated- 
oscilloscope  measurement,  the  video  output  voltage  (proportional  to  log 
power)  could  be  measured.  The  injected  RF  calibration  pulse  would  be 
delayed  to  the  sample  gate  position  desired.  The  RF  signal  source  then 
could  be  step  attenuated  (monitored  through  a  directional  coupler  and 
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calibrated  RF  power  meter)  starting  from  a  level  driving  the  receiver 
into  saturation.  Then  the  "linearity"  of  the  log  IF  response  could  be 
checked  together  with  the  fixed  gain  receiver  dynamic  range  by  stepped 
attenuation  until  the  noise  level  (MDS)  is  reached.  Adjustments  on  the 
delay  of  the  injected  pulse  should  reproduce  the  slope  of  the  STC  cir¬ 
cuitry  (dB/us)  built  into  the  receiver  in  the  fixed  STC  "weather"  mode 
of  operation. 

2.5.4  Data  Recording  System 

The  log  amplitude-constant  gain  characteristic  of  the  video  channel 
is  favorable  for  on-board  gating,  processing,  and  recording  of  the  back- 
scattered  volumetric  returns.  The  video  voltage  would  undergo  range 
gating,  boxcarring  (sample  and  hold),  A/D  conversion,  buffer  storage,  and 
integration  followed  by  integrator  dumping,  formatting  and  digital 
recording.  Timing  circuitry  will  be  needed  to  keep  track  of  integration 
counts,  to  provide  time  marks,  and  to  identify  the  records  on  tape.  The 
level  of  transmitter  power  (probably  also  variable  with  time)  will  also 
need  to  undergo  similar  processing  and  recording;  however,  the  trans¬ 
mitted  power  measurement  would  require  an  A/D  converter  with  a  much 
smaller  dynamic  range.  The  exact  configuration  and  specification  of  the 
equipment  will  depend  on  the  outcome  of  preliminary  tests.  For  example, 
if  the  receiver  dynamic  range  is  determined  to  be  100  dB,  the  size  of 
the  A/D  converter  required  for  better  than  \  dB  quantization  of  amplitude 
could  be  8  bits  (i.e.,  2n  -  1  =  255  levels  at  0.392  dB/level).  Care 
should  also  be  taken  to  provide  for  system  calibration  recordings  at 
frequent  intervals  prior  to,  during,  and  after  mission  flights  (depending 
on  system  stability)  to  assure  meaningful  and  repeatable  reflectivity 
measurements . 

2.5.5  Conclusions  and  Recommendations 

The  use  of  an  on-board  weather  radar  to  collect  airborne  measurements 
of  liquid  water  content  appears  to  be  feasible  providing  that  equipment 
sensitivity  and  stability  are  proven  to  be  adequate.  The  narrow  pencil 
beam  and  high  range  resolution  capability  of  the  AN/APQ-122  (V)  8  systems 
operating  at  Ka-band  (%33  GHz)  and  X-band  (-v9.4  GHz)  is  particularly 
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suited  to  these  requirements.  The  versatility  of  flying  at  designated 
heights  and  locations  is  desirable  as  is  the  availability  for  direct 
correlation  with  other  on-board  sensor  data.  The  sensing  of  liquid  water 
content,  though  indirect,  is  derivable  through  radar  reflectivity  back- 
scatter  measurements  in  a  volume  ahead  of  the  aircraft.  The  accuracy  to 
which  this  can  be  accomplished  depends  to  a  large  part  on  the  care  with 
which  the  systems  are  calibrated  and  also  on  the  equipment  design  and 
stability. 

It  is  believed  that,  through  the  recommended  series  of  system  checks 
enough  of  the  system  characteristics  may  be  determined  to  provide  input 
to  the  design  of  interface  equipment  for  an  on-board  digital  recording 
system  of  pre-integrated  data.  It  does  not  appear  that  any  modification 
of  the  radar  system  will  be  required  if  the  outcome  of  the  preliminary 
tests  proves  favorable.  The  design  of  the  interface  equipment  and  record¬ 
ing  system  could  proceed  following  an  analysis  of  the  preliminary  data 
and  the  availability  of  the  performance  specifications. 

2.6  Task  6  -  Cloud  and  Precipitation  Descriptions  for  Ten 

Case  Studies 

The  synoptic  and  mesoscale  cloud  and  precipitation  patterns  for  nine 
cases  at  Wallops  Island,  Virginia  and  one  case  at  Kwajalein  were  prepared 
Cloud  and  precipitation  cross  sections  were  drawn  using  all  the  synoptic, 
satellite  and  special  data  which  were  available  in  the  specific  region. 
These  analyses  were  prepared  for  AFGL  scientists  who  used  the  data  as 
part  of  their  continuing  research  on  the  growth  of  precipitation. 
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PROGRAM  CLOUD7T 


1.  DESCRIPTION 

Program  CLOUD7T  unpacks  the  ALCOR  data  tape  (ADT)  to  creat  a 
7-track  tape  for  the  AFGL  weather  program.  It  can  also  be  set  up  to 
generate  a  pulse  listing  giving  the  time  to  milliseconds  of  any 
pulse  contained  in  the  particular  ALCOR  data  tape. 

CLOUD7T  produces  one  output  tape  known  as  the  "weather"  tape 
which  can  be  processed  in  terms  of  the  amplitude  or  the  phase.  It 
also  can  print  a  listing  of  pulse  numbers  along  with  their 
corresponding  times. 

The  input  to  CL0UD7T  consists  of  a  9-track  ALCOR  data  tape. 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS  (CL0UD7T  in  UPDATE  FILE  under  TUNG  ID) 

2.1.1  TO  LIST  PULSE  INFORMATION 

JOBNM,CM65000,T300,NT1 , STMFB .  ACT  NO.  NAME 

VSN, TAPE1=MXXXXX/HD. 

ATTACH , OLDPL , CLOUD7T , TD=TUNG . 

UPDATE, N. 

FTN, I=COMPILE,SL,PL=7777777. 

REQUEST, TAPE 1, HD, L, NT, NR.  { MXXXXX/NORING) 

LDSET (PRESET=ZERO) 

LGO. 

7/8/9 

* IDENT  CHECK 
♦DELETE  CLOUD7T . 39 
L IST=1 
*C  CL0UD7T 
7/8/9 

DATA  CARDS 
fi/7/8/9 

(In  the  program  CLOUD7T ,  line  39  reads  LIST=0;  to 
pulse  list  it  must  read  LIST=1.  The  program  is  set  up  to 
tenth  pulse.  This  may  be  altered  by  changing  the  counter 
339  which  is  set  at  10  but  may  be  any  number.) 

2.1.2  TO  CREAT  7-TRACK  AMPLITUDE  TAPE 

JOBNM,CM65000,T500,NT1,MT1, STMFB.  ACT  NO.  NAME 
VSN, TAPE 1=MXXXXX/HD. 

VSN , TAPE 1 5=LYCXXX . 

ATTACH, OLDPL, CL0UD7T, ID=TUNG. 

UPDATE (F) 

FTN, I =COMPI LE ,SL,PL=7777777. 

1 

1 


obtain  a 
list  each 
in  line 


REQUEST, TAPE1, HD, L,NT,NR.  ( MXXXXX/NORING) 

REQUEST, TAPE15, HY, RING.  (LYCXXX/RING) 

LDSET (PRES£T=ZERO) 

LGO. 

7/8/9 

*C  CLOUD7T 
7/8/9 

DATA  CARDS 
6/7 /8/9 

(No  alterations  in  the  CL0UD7T  program  are  needed) 
2.1.3  TO  CREAT  7-TRACK  PHASE  TAPE 

JOBNM , CM 6 50 00 , T900 , NT 1 , MT1 , STMFB .  ACT  NO.  NAME 
VSN, TAPE1=MXXXXX/HD. 

VSN, TAPE1 5=LYCXXX. 

ATTACH , OLDPL , CLOUD7T , TD=TUNG . 

UPDATE, N. 

FTN , I=COMPTLE, SL, PL=7 777777. 

REQUEST, TAPE  1 , HD, L, NT, NR.  (MXXXXX/NORING) 

REQUEST , TA  PE  1 5 , HY, RING.  (LYCXXX/RING) 

LDSET (PRESET=ZFRO) 

LGC. 

7/P/9 

*  IDENT  CHECK 
♦INSERT  CLOUD7T.6 

DIMENSION  TG  (340) 

♦INSERT  CLOUD7T.1? 

DIMENSION  XPTBL (255) ,XLPH(170) ,XRPH(170) 
♦INSERT  CLOUD7T.79 
N  =  0 

DO  21  K=1 , 255 
N=N+ 1 

21  XPTBL(N)=XNBUF (K) 

♦INSERT  CLOUD7T. 301 
DO  38  K  =  1 ,  170 
J=K  +  170 
N=I LCPHA (K )  +  128 
XL  PH  (K )  =XPTB  L  (N  ) 

NCON (K ) =XLPH (K) *100. 

M=TRCPHA (K)  +  128 
XRPH (K ) =XPTB  L ( M) 

NCON  (J )  =XRPH  (K )  *  1 00  . 

38  CONTINUE 
♦INSERT  CLOUD7T .325 

CALL  GBYTES (NOUT(l)  ,NCON(l)  ,1,14,1,170) 

CALL  GBYTES (NOUT(l)  , TG  ( 1  )  ,  0 ,  1  ,  1 4 , 340 ) 

DO  80  1=1,340 
I F ( TG ( I ) . EQ. 0 )  GO  TO  80 
TG  ( I )  =  . NOT .  NCON ( I ) 

CALL  GBYTE  ( IG  ( I )  , NCON ( I )  ,  46 , 1 4 ) 

NCON  (!)=(-]  )  *NCON  ( I ) 
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80  CONTINUE 

*  INSERT  CLOUD7T . 480 

DIMENSION  ISG (170) 

*  INSERT  CL0UD7T. 503 

NSKT  P=NSO  +  8*171 
CALL  SKIP (NO,N,NSKIP) 

CALL  GBYTES (IN (N) , ISG(l) , NSKI P, 1 , 7 , 1 7 0 ) 
NSKI P=NSO  +  8*171  +  1 
CALL  SKIP(NO,N, NSKI P) 

CALL  GBYTES (IN (N ) ,IDAT(171) ,NSKIP, 7, 1 , 170) 
DO  23  IG=171 , 340 
NG=IG-1 70 

2  3  I  DAT ( IG ) = (-1 ) **TSG (NG ) *  I  DAT ( IG ) 

♦INSERT  CLOUD7T . 509 

NSKI P=NSO  +  8*522 
CALL  SKIP(NO,N,NSKIP) 

CALL  GBYTES (IN (N) ,  TSC(l) , NSKI P , 1 , 7 , 1 7 0 ) 
NSKI P=NSO  +  8*522  +  1 
CALL  SKIP (NO,N,NSKIP) 

CALL  GBYTES (IN (N) , I DAT (5 1 1 ) , NSKI P, 7 , 1 , 1 70 ) 
DO  24  IG=5 11,580 
NG  =  IG-2  50 

24  I  DAT ( IG ) = (-1)**ISC(NG)*I DAT ( IG ) 

*C  CLOUD7T 
7/8/9 

DATA  CARDS 
6/7 /8/9 


2.2 

INPUT 

DATA  CARDS 

(for  all 

CLOUD7T  jobs) 

Card 

#1  (2110, 2F10. 3, 10A4) 

cc 

1-10 

(7200724) 

ILNCH 

Liftoff (Millis) 

cc 

1  1-20 

(1) 

NVALS 

Number  of  start-stops 

cc 

33-40 

(14.79) 

GTESP 

Inter  Gate  Spacing 

in  meters 

cc 

41-80 

ANT-3A 

RUN  ID 

Up  to  40  character 

title 

Ca  rd 

02  (8110) 

cc 

1-10 

(22401) 

NSTART 

Start  Pulse  #  (right  justified) 

cc 

11-20 

(42400) 

NSTOP 

End  Pulse  #  (right 

justified) 

3.  TIME  AND  MEMORY  REQUIREMENTS 


The  65k  in  octal  memory  is  required  to  run  the  program.  The 
500  cp  seconds  is  enough  to  creat  the  amplitude  weather  tape  for 
10000  pulses.  It  requires  900  cp  seconds  to  creat  the  phase  weather 
tape  for  10000  pulses. 

4.  TAPES 

TAPE  1  is  an  INPUT  tape  consisting  of  a  9-track  ALCOR  data  tape 
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from  Lincoln  Laboratory.  The  FORMAT  of  this  tape  is  given  in  ALCOR 
DATA  USERS  AND  RADAR  OPERATORS  MANUAL. 

TAPE15  is  an  OUTPUT  tape  consisting  of  a  7-track  800  BPI 
weather  tape.  The  FORMAT  of  this  tape  is  listed  in  TABEL  A-l.  A 
maximum  of  10000  pulses  (records)  can  be  stored  on  the  regular 
7- track  tape. 

5. SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

REFC ( E , R , DEE , DRP) 

READGS ( I  EOF , I  ERR) 

HEADT ( ISIG, IEQM) 

CALIBR (XNBUF) 

UNPACK (LIST) 

SKIP  (NO) 

PPAGC (NO) 

OPAGC (NO) 

REAL  360  (N  ,  B) 

CHAR  36 0  ( INT , N) 

PACKK 

GBYTES 
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PROGRAM  PRHISEC 


1.  DESCRIPTION 

Program  PRHISEC  processes  the  weather  tape  created  from 
CL0UD7T  to  produce  the  line  printer  RHI  plots  of  reentry  weather 
scans.  The  scale  of  the  plot  is  controlled  by  the  data  cards  and  it 
is  possible  to  obtain  several  different  plots  with  different 
altitude  vs.  range  of  the  same  data  at  the  same  time. 

Input  to  the  PRHISEC  program  includes  one  or  more  weather 
tapes  previously  obtained  from  the  CL0UD7T  program. 

The  PRHISEC  program  is  set  up  to  use  a  tracking  gate  of  52.5 
(line  73)  and  a  correlation  factor  of  82.5  (line  81).  These  values 
may  change  for  different  missions. 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS  (PRHISEC  program  in  UPDATE  FILE  under  TUNG  ID) 

JOBNM ,T700,CM]70000,MT2.  ACT  NO.  NAME 

VSN, TAPE! =LYCXXX . 

VSN, TAP£2=LYCXXX . 

ATTACH, OLDPL, PRHISEC, ID=TUNG , MR=] . 

UPDATE (F) 

FTN , I=COMPILE. 

REQUEST, TAPE1 ,HY, L.  ( LYCXXX/NOR ING ) 

REQUEST, TAPE2, HY, L.  ( LYCXXX/NOR ING ) 

LDSET ( PRESET=Z ERO ) 

LGO. 

REWIND, TAPES . 

C0PYCF,TAPE6, OUTPUT. 

REWIND, TAPES . 

COPYCF , TAPES , OUTPUT. 

7/8/9 

*C  PRHISEC 
7/8/9 

DATA  CARDS 
S/7/8/9 

(This  job  is  set  up  for  the  case  where  two  weather  tapes  are 
needed  as  data  input.  If  only  one  weather  tape  is  needed  only  TAPE  1 
would  be  used  with  the  VSN  and  REQUEST  cards  and  the  number  of 
tapes  would  *">e  changed  on  the  JOB  card.  Currently  PRHISEC  is  set  up 
to  handle  no  more  than  two  weather  tapes.  Notice  also  that  the 
output  tape  (TAPES)  is  rewound  and  copied  twice.  This  will  result 
in  two  sets  of  plots  being  generated.  If  only  one  set  is  desired 
remove  one  REWIND  and  COPYCF  of  TAPES.) 


2.2  INPUT  DATA  CARDS 


The  input  data  cards  are  in  NAMELIST  form  (SPARAM) .  For  normal 
processing,  one  plot  is  generated  with  the  desired  Height  vs. 
Range.  The  data  cards  are  set  up  as  follows: 

SPARAM 

DBMIN=0 .,XMIN=42.0,XMAX=50.0,YMIN=2.00,YMAX=8.00, 

BTIME=0 . ,  NV=170 . ,  ETIME=0 . , 

LASPLS=193768 , 

N0NPLS=0 , 

Send 

SPARAM 

DBMIN=999. 9, 

Send 


The  second  set  of  NAMELIST  is  set  to  terminate  the  job.  If 
desired,  more  plots  with  different  Height  and  Range  values  can  be 
generated  by  adding  more  sets  of  NAMELIST  PARAM. 


2. 3  NAMELIST  PARAM 

VARIABLES 

NAME 

TYPE 

MEANING 

DBM  IN 

R 

Current  not  used,  set  to  1.0. 

DBMIN  set  to  999.9  for  end  of  job. 

XM  IN 

R 

Minimum  horizontal  range  of  plot  in  km. 

XMAX 

R 

Maximum  horizontal  range  of  plot  in  km. 

YM  IN 

R 

Minimum  vertical  height  of  plot  in  km. 

YMAX 

R 

Maximum  vertical  height  of  plot  in  km. 

BTIME 

R 

Time, in  total  seconds  of  day,  at  which 
to  begin  processing  of  scan. 

ETIME 

R 

Time  to  end  processing  of  scan,  to  process 
whole  scan,  set  BTIME  and  ETIME  to  0. 

NV 

I 

Number  of  data  values  per  tape  record, 
current  set  to  ]70. 

LASPLS 

I 

Last  pulse  number  of  scan,  to  end 
processing  of  data. 

NONPLS 

I 

Pulse  number  of  preceeding  scan  on  tape, 
set  to  0  if  scan  starts  at  beginning  of  tape 

3.  TIME  AND  MEMORY  REQUIREMENT 


The  117k  in  octal  of  memory  is  required  to  run  this  program. 
Approximately  300  cp  seconds  processing  time  is  required  for  one 
scan  plot. 

«t.  TAPES 

One  or  two  amplitude  weather  tapes  (format  in  TABLE  A-l) 
previously  processed  by  CL0UD7T  are  needed  as  input  to  run  this 
program.  TAPE1  and  TAPE2  are  reserved  for  the  input  tapes.  TAPE5  is 
for  the  NAMELIST  input  data  cards.  TAPES  is  for  the  printed 
output, the  deck  is  set  up  to  rewind  and  copy  to  the  output  file. 


A-S 


5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

READK ( I  END) 

GBYTE 

GBYTES 
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PROGRAM  TRADEX 


1.  DESCRIPTION 

Program  TRADEX  was  developed  to  process  the  AFGL  SPECIAL 
TRADEX  weather  tape  to  generate  line  printer  plots  for  RHI  maps. 

The  input  weather  tape  was  specially  created  by  Lincoln 
Laboratory  for  AFGL  use.  The  header  record  of  the  tape  provides  the 
basic  information  for  processing.  The  output  plot  contains  a  map  of 
DBZ . 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS 

JOBNM,T300,CM120000, TP2 .  ACT  NO.  NAME 
FTN, SL, PL=500000. 

VSN , TAPE  1 =TR0880 . 

VSN, TAPE2=TR0940 . 

REQUEST, TAPE1 , HY , NR .  (TR0880 /NORING) 

REQUEST, TAPE2,HY, NR.  (TR0940/NORING) 

FILE (TAPE  1 , BT= I , RT=W,MRL=51 30 ,MBL=51 30 ) 

FILE (TAPE2,BT=I , RT=W , MRL=5 1 30 , MBL=5 1 30 ) 

LDS  ET (FI LES=TAPE 1 ) 

LDSET ( FILES=TAPE2) 

LDSET (PRESET=ZERO) 

LGO . 

REWIND, TAPES. 

COPYCF , TAPES , OUTPUT . 

7/8/9 

SOURCE  DECK 
7/8/9 

DATA  CARDS 
S/7/8/9 

2.2  TNPUT  DATA  CARDS 

The  input  data  cards  are  in  NAMELIST  form  ($PARAM) .  For  normal 
processing,  one  plot  is  generated  with  the  desired  height  vs. 
range.  The  data  cards  are  set  up  as  follows: 

SPARAM 

DBM  TN= 1 . 0,XMIN=50.0,XMAX=80.0,YMIN=0.0,YMAX=20.0, 

BTIMF.=0 . 0  ,  ETIM  E=0  .  0 ,  NV=500, 

SEND 

SPARAM 

DBM IN=999 . 9 , 

SEND 

The  second  set  of  NAMELIST  is  set  to  terminate  the  job.  If 
desired,  more  plots  with  different  height  and  range  scales  can  be 
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generated  by  adding  more  sets  of  NAMELIST  PARAM. 
2.3  NAMELIST  PARAM  VARIABLES 


NAME 

TYPE 

MEANING 

DBM  IN 

R 

Current  not  used,  set  to  1.0. 

DBM IN  set  to  999.9  for  end  of  job. 

XM  IN 

R 

Minimum  horizontal  range  of  plot  in  km. 

XMAX 

R 

Maximum  horizontal  range  of  plot  in  km. 

YM  IN 

R 

Minimum  vertical  height  of  plot  in  km. 

YMAX 

R 

Maximum  vertical  height  of  plot  in  km. 

BTIME 

R 

Start  time,  set  to  0.0. 

The  header  record  of  input  tape  will  provide 
the  tape  start  time  in  GMT  total  seconds. 

ETIME 

R 

Stop  time,  set  to  0.0. 

The  header  record  of  input  tape  will  provide 
the  tape  stop  time  in  GMT  total  seconds. 

NV 

I 

Number  of  range  gates  per  pulse  (number 
of  data  values  per  tape  record).  Set 
to  500  for  current  use.  500  is  the 
maximum  array  dimension  in  the  program. 

The  header  record  of  the  input  tape 
will  provide  the  value  for  each  mission. 

3.  TIME  AND  MEMORY  REQUIREMENT 


The  120k  in  octal  of  memory  is  required  to  run  this  program. 
Approximately  300  cp  seconds  processing  time  is  required  for  one 
scan  plot. 

A.  TAPES 

TAPE!  and  TAPE 2  are  reserved  as  the  input  tapes  called  AFGL 
SPECIAL  TRADEX  weather  tape  (format  in  TABLE  A-2)  from  Lincoln 
Laboratory.  TAPE5  is  for  the  NAMELIST  input  data.  TAPES  is  for  the 
printed  output,  the  deck  is  set  up  to  rewind  TAPES  and  copy  to  the 
output  file. 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

READT ( I  END) 

GBYTE 

GBYTES 
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PROGRAM  APLOT  AND  TPLOT 


1.  DESCRIPTION 

Programs  APLOT  and  TPLOT  were  developed  to  process  the 
vertical  (zenith)  scan  radar  data.  They  creat  a  plot  from  CALCOMP 
plotter  to  plot  DB7.  values  vs.  height  in  kilometers.  Program  APLOT 
processes  the  ALCOR  radar  weather  tape  created  by  the  CL0UD7T 
program  and  program  TPLOT  processes  the  AFGL  Special  TRADEX  weather 
tape  provided  by  Lincoln  Laboratory. 

Input  consists  of  the  ALCOR  or  TRADEX  weather  tape.  The  output 
is  an  eleven  inches  by  eleven  inches  online  plot  and  a  listing  of 
height  with  DBZ  values. 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS 

JOBNM , T700 , MT2 .  ACT  NO.  NAME 

FTN ,SL,PL= 50000. 

VSN (TAPE  3=LYCXXX ) 

VSN (TAPE4=LYCXXX) 

REQUEST, TAPE3,HY, NR.  ( LYCXXX/NOR ING ) 

REQUEST, TAPE4,HY, NR.  (LYCXXX/NOR ING) 

COPYBF,TAPE3,TAPEl . 

COPYBF,TAPE4,TAPEl. 

REWIND, TAPE1. 

REQUEST, PLOT, *Q. 

ATTACH, PEN, ONLINEPEN. 

LIBRARY, PEN. 

DISPOSE, PLOT, *PL. 

LDSET, FILES=TAPE1 , PRESET=ZERO. 

LGO. 

7/8/9 

SOURCE  DECK 
7/8/9 

DATA  CARDS 
S/7/8/9 

NOTE:  Program  TPLOT  needs  to  insert  a  FILE  card  in  contral  cards. 
FILE (TAPE  1 , BT=T , RT=W,MRL=51 30, MB L* 5130) 

2.2  INPUT  DATA  CARDS 

The  input  data  cards  are  in  NAMELIST  form  ($PARAM) .  They  are 
set  up  as  follows: 

$PARAM 

YMTN=5.0,YMAX=20.0,BTIME=0. ,ETIME=0. ,NV=170, 

LASPLS=XXXXXX , 

$END 


SPARAM 

DBMIN=999 . 9 , 
SEND 


2.2.1  NAMELIST  PARAM  VARIABLES 


NAME 

TYPE 

MEANING 

YM  IN 

R 

Minimum  vertical  range  of 

plot  ( km) . 

YMAX 

R 

Maximum  vertical  range  of 

plot  ( km) . 

BTIME 

R 

Beginning  time. 

ETIME 

R 

End  time. 

NV 

I 

Number  of  data  values  per  tape  record, 

170  for  APLOT.  Value  changed  of  different 
mission  for  TPLOT,  500  is  the  maximum 
array  dimension. 

LASPLS 

I 

Last  pulse  number  of  data 
be  processed . 

to 

DBMIN 

R 

Set  to  999.9  for  end  of  processing. 

3.  TIME  AND  MEMORY  REQUIREMENT 

This  program  can  be  run  under  the  default  core.  Approximately 
S00  cp  seconds  processing  time  is  required  to  process  an  input  tape 
containing  10000  pulses. 

4.  TAPES 

TAPE1  is  an  input  file  copied  from  the  input  tapes  (TAPE3  and 
TAPE* ) .  For  program  APLOT,  the  input  tapes  are  the  ALCOR  amplitude 
weather  tapes  previosly  created  from  the  CL0UD7T  program  (format  in 
TABLE  A-l ) .  For  program  TPLOT,  the  input  tapes  are  the  special  AFGL 
TRADEX  weather  tapes  (format  in  TABLE  A-2)  provided  by  Lincoln 
Laboratory.  TAPE5  is  used  for  the  NAMELIST  input.  TAPES  is  used  for 
the  printer  output. 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS 


5.1  APLOT 

READK (IEND) 

GBYTE 

GBYTES 

5.2  TPLOT 

READT ( I  END) 

GBYTE 

GBYTES 
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PROGRAM  MORT2 


1.  DESCRIPTION 

The  MORT2  program  was  modified  from  the  PRHISEC  program  to  get 
one  second  average  information.  The  printout  of  this  program  is 
similar  to  the  MOIST  printout  from  Lincoln  Laboratory.  It  contains 
time,  range  (km),  altitude  (km),  azimuth  (radians)  ,  elevation 
angle  (radians) ,  reflectivity  (DBZ)  and  radar  cross  section  for 
each  second.  If  desired  it  can  also  be  set  up  to  creat  an  output 
tape  in  MOIST  TAPE  format. 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS  (M0RT2  in  UPDATE  FILE  under  TUNG  ID) 

JOBNM,T300,CM1200C0,TP2.  ACT  NO.  NAME 

VSN , TAPE  1 1 =LYCXXX . 

VSN, TAPE! 2=LYC XXX. 

REQUEST, TAPE  1 1 , HY, NR.  ( LYCXXX/NORING) 

REQUEST, TAPE  12, HY, NR.  (LYCXXX/NORING) 

ATTACH, 0LDPL,M0RT2, ID=TUNG . 

UPDATE (F) 

FTN , I=COMPILE,SL,PL=7777777. 

LDSET (PRESET=ZERO) 

LGO. 

7/8/9 
*C  MORT2 
7/8/9 

DATA  CARDS 
S/7/8/9 

(This  job  is  set  up  for  the  case  where  two  weather  tapes  are 
needed  but  as  many  as  ten  tapes  may  be  used  as  input  data  if 
needed.  Whatever  number  input  tapes  are  used,  the  first  tape  should 
be  assigned  TAPE11;  the  second  TAPE12;  etc.,  up  to  TAPE21,  being 
certain  to  include  a  set  of  VSN  and  REQUEST  cards  for  each  and 
denoting  the  total  number  on  the  JOB  card.) 

2.2  INPUT  DATA  CARDS 

The  input  data  cards  are  in  NAMELIST  form  (SPARAM).  The  data 
cards  are  set  up  as  follows: 

SPARAM 
NT A PE  =  2, 

DBM IN= 1.0, XM IN=0 . 0 , XMAX=90 . 0 , YM TN=0 . 0 , YMAX=2 0 . 0 , 

BTTM  E=0 . ,  NV=1 70 ,  ETIME=0 . , 

NONPLS=0 , 

LASPLS  =  382  32 , 

SEND 

SPARAM 

DBM IN=999 . 9 , 
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SEND 


(The  second  set  of  NAMELIST  is  set  to  terminate  the  job.) 


2.3  NAMELIST 

PARAM  VARIABLES 

NAME 

TYPE 

MEANING 

NTAPE 

I 

Number  of  input  tapes  used  in  job. 

DBM  IN 

R 

Current  not  used,  set  to  1.0. 

Set  to  999.9  for  end  of  job. 

XM  IN 

R 

Minimum  horizontal  range  of  data  (km). 

XMAX 

R 

Maximum  horizontal  range  of  data  (km). 

YM  IN 

R 

Minimum  vertical  height  of  data  (km). 

YMAX 

R 

Maximum  vertical  height  of  data  (km). 

BTIME 

R 

Time  (total  seconds  of  day)  at  which 
to  begin  processing. 

ETIME 

R 

Time  (total  seconds  of  day)  at  which 
to  end  processing.  To  process  all  data 
set  BTIME  and  ETIME  to  0. 

NV 

I 

Number  of  data  values  per  tape  record, 
current  is  set  to  170. 

NONPLS 

I 

Pulse  number  of  preceeding  data  on 
tape  (set  to  0  if  processing  starts 
at  beginning  of  tape). 

LASPLS 

I 

Last  pulse  number  of  processing, 
to  end  processing  of  data. 

3.  TIME  AND  MEMORY  REQUIREMENT 

The  120k  in  octal  of  memory  is  required  to  run  this  program. 
Approximately  60  cp  seconds  processing  time  is  required  to  run  the 
data  on  one  weather  tape. 

4.  TAPES 

TAPE 1 1  to  TAPE 21  are  reserved  for  the  input  tapes  (format  in 
TABLE  A- 1 )  .  TAPE3  is  an  output  tape  in  MOIST  TAPE  format. 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

READK  ( I  END) 

GBYTE 

GBYTES 
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PROGRAM  PPOP 


1.  DESCRIPTION 

Program  PPOP  processes  the  amplitude  weather  tape  to  produce 
time  series  maps  for  each  second  of  time.  The  program  creates  four 
output  sets  consisting  of:  the  mapping  for  the  amplitude  of 
principal  polarization  (PP);  the  mapping  for  the  amplitude  of 
orthogonal  polarization  (OP);  the  mapping  for  the  difference 
between  the  principal  and  orthogonal  polarization  (PP-OP) ;  and  the 
listing  for  the  height  (km)  of  the  first  and  last  range  gate  used. 
PPOP  uses  a  six-gate  averaging  scheme  with  the  last  (29th)  value 
being  a  two-gate  average.  Limitations  in  page  size  (OUTPUT)  with 
intent  on  obtaining  the  smallest  number  of  gates  to  be  averaged  led 
to  this  averaging  scheme.  The  values  of  principal  and  orthogonal 
polarization  and  their  difference  are  mapped  in  tenths  DBZ  with 
their  appropriate  sign  (blank  equals  positive). 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS  (PPOP  in  UPDATE  FILE  under  TUNG  ID) 

JOBNM ,  T500 ,  CM  1 20000 , TP1 .  ACT  NO.  NAME 

VSN, TAPE 1 =LYCXXX . 

ATTACH, OLDPL, PPOP, ID=TUNG . 

UPDATE (F) 

FTN , I =COMPI LE ,SL,PL=7777777. 

REQUEST, TAPE 1,HY, NR.  ( LYCXXX/NOR ING ) 

LDSET (PRESET=ZERO) 

LGO. 

REWIND, TAPE*. 

COPY, TAPES, OUTPUT. 

REWIND, TAPE7. 

COPY, TAPE7, OUTPUT. 

REWIND, TAPER. 

COPY, TAPE8, OUTPUT. 

REWIND, TAPE3. 

COPY, TAPE  7, OUTPUT. 

7/8/9 
*C  PPOP 
7/8/9 

DATA  CARDS 
S/7/8/9 

(If  a  second  tape  is  needed  as  input  it  would  be  designated 
TAPE2  and  a  VSN  and  REQUEST  card  for  TAPE2  would  be  added  following 
the  present.  VSN  and  REQUEST  cards.  Also,  the  number  of  tapes 
indicated  on  the  JOB  card  would  be  changed  to  2.) 

2.2  INPUT  DATA  CARDS 

The  input  data  cards  are  in  NAMELIST  form  ($PARAM)  and  they 
are  set  up  as  follows: 
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SPARAM 

DBMIN=1.0,XMIN=0.0,XMAX=90.0,YMIN=0.0,YMAX=20.0, 
BTIME=0 . ,  NV=170,  ETTME=0 . , 

I  PASS=7 , 

NONPLS=0, 

LASPLS=44800, 

SEND 

SPARAM 

DBM IN=999 . 9 , 

SEND 

2.7  NAMELIST  PARAM  VARIABLES 


NAME 

TYPE 

MEANING 

DBM  IN 

R 

Current  not  used,  set  to  1.0. 

DBMIN  set  to  999.9  for  end  of  job. 

XM  IN 

R 

Minimum  horizontal  range  of  data  (km). 

XMAX 

R 

Maximum  horizontal  range  of  data  (km). 

YM  IN 

R 

Minimum  vertical  height  of  data  (km). 

YMAX 

R 

Maximum  vertical  height  of  data  (km). 

BTIME 

R 

Time,  in  total  seconds  of  day,  at  which 
to  begin  processing  scan. 

ETIME 

R 

Time  to  end  processing  of  data,  to 
process  all  of  data  set  BTIME  and 

ETIME  to  0. 

NV 

I 

Number  of  data  values  per  tape  record, 
current  set  at  170. 

I  PASS 

I 

Identification  number  of  pass  working  on 

NONPLS 

I 

Pulse  number  of  preceeding  data  on  tape, 
set  to  0  if  data  starts  at  beginning  of 

LASPLS 

I 

Last  pulse  number  of  data,  to  end 
processing  of  data. 

3.  TIME  AND  MEMORY  REQUIREMENT 


The  120k  in  octal  of  memory  is  required  to  run  this  program. 
Approximately  500  cp  seconds  processing  time  is  needed  to  process 
two  input  tapes  (each  about  10000  pulses)  of  data. 

4.  TAPES 

TAPE 1  and  TAPE 2  are  reserved  for  input  tapes  (format  in  TABLE 

A-l )  . 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

READK  (IEND) 

GBYTE 

GBYTES 
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PROGRAM  PPOPF 


1.  DESCRIPTION 

The  PPOPF  program  processes  the  phase  weather  tape  to  produce 
an  every  second  phase  difference  listing.  It  lists  the  average 
phase  difference  in  radians  for  all  170  gates  along  with  tracking 
gate  (52)  information.  The  tracking  gate  information  consists  of: 
the  time  to  milliseconds;  the  true  range  (km);  the  azimuth 
(radians);  the  elevation  angle  (radians);  the  altitude  (km);  the 
range  with  respect  to  the  X-axis  (km);  the  range  with  respect  to 
the  Y-axis  (km);  and  the  phase  difference.  The  170  values  listed 
are  the  one  second  phase  differences  of  the  170  range  gates. 

Input  to  the  PPOPF  program  consists  of  one  or  two  phase 
weather  tapes  previously  processed  by  the  CLOUD7T  program. 

Output  from  the  PPOPF  program  consists  of  a  listing  of  the 
phase  difference  for  all  170  gates  of  each  200  pulse  average  or, 
one  second  of  time. 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS  (PPOPF  in  UPDATE  FILE  under  TUNG  ID) 

JOBNM ,T500,CM120000,TP1 .  ACT  NO.  NAME 
VSN, TAPE 1=LYC XXX. 

ATTACH , OLDPL , PPOPF , I D=TUNG . 

UPDATE (F) 

FTN, I=COMPILE,SL,PL=7777777. 

REQUEST, TAPE] , HY , NR .  (LYCXXX/NORING) 

LDSET ( PRESET=Z ERO ) 

LGO. 

7/8/9 
*C  PPOPF 
7/8/9 

DATA  CARDS 
5/7/R/9 

(If  a  second  tape  is  needed  as  input  it  would  be  designated 
TAPE2  and  a  VSN  and  REQUEST  card  for  TAPE2  would  be  added  following 
the  present  VSN  and  REQUEST  cards.  Also,  the  number  of  tapes 
indicated  on  the  JOB  card  would  be  changed  to  2.) 

2.2  TNPUT  DATA  CARDS 

The  input  data  cards  are  in  NAMELIST  form  (SPARAM)  and  they 
are  set  up  as  follows: 

SPARAM 

DBM IN= 1 . 0,XMTN=0.0, XMAX=90 . 0 , YM IN=0 . 0 , YMAX=20 . 0 , 

BTIME=0 „ ,  NV= 170,  ETIME=0 . , 

NONPLS  =  0 , 
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LASPLS  =  4  4  800 , 

SEND 

SPARAM 

DBMIN=999. 9, 

SEND 

2.2  NAMELIST  PARAM  VARIABLES 


NAME 

TYPE 

MEANING 

DBM  IN 

R 

Current  not  used,  set  to  1.0. 

DBMIN  set  to  999.9  for  end  of  job. 

XM  IN 

R 

Minimum  horizontal  range  of  data  (km). 

XMAX 

R 

Maximum  horizontal  range  of  data  (km). 

YM  IN 

R 

Minimum  vertical  height  of  data  (km). 

YMAX 

R 

Maximum  vertical  height  of  data  (km). 

BTIME 

R 

Time,  in  total  seconds  of  day,  at 
which  to  begin  processing  scan. 

ETIME 

R 

Time  to  end  processing  of  data,  to 
process  all  of  data  set  BTIME  and  ETIME 
to  0 . 

NV 

I 

Number  of  data  values  per  tape  record, 
current  set  at  170. 

NONPLS 

I 

Pulse  number  of  preceeding  data  on  tape 
set  to  0  if  data  starts  at  beginning  of 

LASPLS 

I 

Last  pulse  number  of  data,  to  end 
processing  of  data. 

?.  TIME  AND  MEMORY  REQUIREMENT 


The  120k  in  octal  of  memory  is  required  to  run  this  program. 
Approximately  500  cp  seconds  processing  time  is  needed  to  process 
two  input  tapes  (each  about  10000  pulses)  of  data. 

4.  TAPES 

TAPE1  and  TAPE2  are  reserved  for  input  tapes  (format  in  TABLE 

A  — 1 )  . 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

READK  ( I  END) 

GBYTE 

GBYTES 
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PROGRAM  PPOPX 


1.  DESCRIPTION 

Program  PPOPX  was  developed  to  improve  the  data  processing 
techniques  by  using  polarization  measurements  of  hydrometeors  to 
process  ALCOR  data  tapes.  Four  final  terms  for  each  range  gate  of 
every  second  were  computed  and  processed  to  generate  the  RHI  maps. 
These  four  terms  are  reflectivity  in  DBZ ,  cancellation  in  DB , 
correlation  in  percent  and  orientation  in  tens  of  degrees. 

The  input  to  program  PPOPX  include  the  amplitude  weather  tape 
and  the  phase  weather  tape  both  created  by  the  CL0UD7T  program.  The 
output  consists  of  the  four  maps  of  the  four  terms  mapped  for  every 
second . 

2.  DECK  SET  UP 

7.1  CONTROL  CARDS  (PPOPX  program  in  UPDATE  FILE  under  TUNG  ID) 

JOBNM,CM140000,T]000,TP2.  ACT  NO.  NAME 
VSN , TAPE  1 =  LYC  XXX . 

VSN , TAPE  1 1 =LYC  XXX . 

REQUEST , TAPE  1 , HY , NR .  ( LYCXXX/NOR ING ) 

REQUEST, TAPE11, HY,NR.  ( LYCXXX/NOR ING) 

ATTACH, OLDPL, PPOPX, ID=TUNG . 

UPDATE (F) 

FTN, I -COMPILE, SL,PL*7777777. 

LDSET (PRESET=ZERO) 

LGO. 

7/8/9 
*C  PPOPX 
7/8/9 

DATA  CARDS 
fi/7/8/9 

(TAPE)  is  the  amplitude  weather  tape  and  TAPE11  is  the  phase 
weather  tape.) 

2.2  INPUT  DATA  CARDS 

The  input  data  cards  are  in  NAMELIST  form  ( SPARAM )  and  the 
data  cards  are  set  up  as  follows: 

SPARAM 

DBM  IN- 1 . 0 , XMTN-S3. 0,XMAX  =  7  3. P, YMIN-7. 5, YMAX=1 2.5, 

BTIM  E=0 .  ,  NV= 170,  ETIME-0 . , 

NONPLS=0 , 

LASPLS  =  33400, 

SEND 

$PARAM 

DBMTN-999. 9, 

SEND 
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The  second  set  of  NAMELIST  is  set  to  terminate  the  job.  Due  to 
the  rather  lengthy  size  of  the  run  it  is  not  recommended  to  add  on 
any  additional  plots. 

2.3  NAMELIST  PARAM  VARIABLES 

NAME  TYPE  MEANING 

DBM  IN  R  Current  not  used,  set  to  1.0. 

DBMIN  set  to  999.9  for  end  of  job. 

XM IN  R  Minimum  horizontal  range  of  plot  in  km. 

XMAX  R  Maximum  horizontal  range  of  plot  in  km. 

YM IN  R  Minimum  vertical  height  of  plot  in  km. 

YMAX  R  Maximum  vertical  height  of  plot  in  km. 

BTIME  R  Time,  in  total  seconds  of  day,  at  which 

to  begin  processing  of  scan. 

ETIME  R  Time,  to  end  processing  of  scan,  to  process 

whole  scan,  set  BTIME  and  ETIME  to  0. 

NV  I  Number  of  data  values  per  tape  record, 

current  set  to  170. 

LASPLS  I  Last  pulse  number  of  scan,  to  end 

processing  of  data. 

NONPLS  I  Pulse  number  of  preceeding  scan  on  tape 

set  to  0  if  scan  starts  at  beginning  of  tape. 

3.  TIME  AND  MEMORY  REQUIREMENTS 

The  140k  in  octal  of  memory  is  required  to  run  this  job. 
Approximately  1000  cp  seconds  processing  time  is  required  for  one 
set  of  plots  involving  about  10000  pulses  on  the  weather  tapes. 

4.  TAPES 

TAPE!  is  reserved  to  input  the  amplitude  weather  tape  and 
TAPE  1 1  is  reserved  to  input  the  phase  weather  tape.  The  format  of 
these  tapes  are  given  in  TABLE  A-l .  TAPE  5  is  used  for  the  NAMELIST 
input.  TAPE  *  is  used  for  the  printer  output  by  rewinding  TAPE  6 
and  copying  to  ouput  file. 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS 

READK ( I  END] 

AMP 
PHASE 
SEC SUM 
PLSET 

RHPLOT ( I JK) 

GBYTE 

GBYTES 
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PROGRAM  WRHISEC  FOR  WALLOPS  ISLAND  DATA 


1.  DESCRIPTION 

Program  WRHISEC  produces  line  printer  plots  of  radar  RHI  scans 
for  Wallops  Island  data.  This  program  has  been  modified  for  years. 
All  the  modifications  are  listed  in  program  comments.  The  current 
version  was  modified  in  july  1977  to  include  vertical  cutoff  for 
rain  attennuation . 

2.  DECK  SET  UP 

2.1  CONTROL  CARDS 

JOBNM , CM  107000, T100 , TP1 .  ACT. NO.  NAME 

FTN, SL, PL=50000. 

VSN, TAPE2=W43292. 

REQUEST, TAPE2, HI, S. 

COPYBF , TAPE 2 , TAPE  1 , 1. 

SKI PF , TAPE  2,7, 17, B. 

COPYBF, TAPE2, TAPE! , 1  . 

REWIND, TAPE1. 

UNLOAD, TAPE2. 

RETURN, TAPE2. 

FILE (TAPE  1 ,BT=K, RT=S ,MBL=1  54  0, MRL=1 54  0, RB=1 , EO=T) 

LDSET, FILES=TAPE1 , PRESET*ZERO. 

LGO. 

REWIND, TAPEf. 

COPYCF , TAPER , OUTPUT . 

REWIND, TAPER. 

COPYCF , TAPER , OUTPUT . 

7/8/9 

SOURCE  DECK 
7/8  /9 

DATA  CARDS 
R/7/8/9 

2.2  INPUT  DATA  CARDS 

The  structure  of  the  input  cards  is  ordered  as  follows: 

Set  1  Title  card 

Set  2  NAMELIST  CORR 

Set  3  N  set  NAMELIST  PARAM  (N  is  number  of  files  to  be 
processed  +  1 ) 

2.2.1  TITLE  CARD  (8A10) 

The  TITLE  CARD  contains  the  data  aquisition  date  and  the  begin 
and  end  time  of  the  first  and  last  files  on  the  tape. 

2.2.2  NAMELIST  CORR  VARIABLES 
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NAME 


TYPE 


MEANING 


R  Glover  constant  used  in  computing  rainfall 

attenuation,  set  to  O.fi  for  runs  without 
rainfall  attenuation- 

R  Used  in  computing  oxygen  attenuation,  set 

to  0.0  for  runs  without  oxygen  attenuation. 

R  Constant  for  7- R  relation. 

R  Constant  for  Z-R  relation. 

R  Vertical  height  (KM)  at  and  above  which 

correction  for  rain  attenuation  stops.  If 
RHI  is  completely  rain  attenuated,  set  RAINHT 
to  YMAX . 


2.2.?  NAMELIST  PARAM  VARIABLES 


NAME 

TYPE 

MEANING 

DBM  IN 

R 

Minimum  DB  of  mapping,  normally  set  to  41.0, 

DBMIN  set  to  999.9  for  end  of  processing. 

XM  IN 

R 

Minimum  horizontal  range  of  plot  (KM). 

XMAX 

R 

Maximum  horizontal  range  of  plot  (KM). 

YM  IN 

R 

Minimum  vertical  range  of  plot  (KM),  normally 
set.  to  0.0. 

YMAX 

R 

Maximum  vertical  range  of  plot  (KM),  normally 
set  to  12.0. 

BTIME 

R 

Beginning  time  of  data  to  be  processed,  in  form 
DDD.HHMMSS  where  DDD  is  julian  day. 

ETIME 

R 

End  time  of  date  to  be  processed.  To  process  all 
data  in  the  file,  set  BTIME  and  ETIME  to  0.0. 

CEL 

R 

Elevation  angle  correction,  set  to  0.0  for  no 
correction . 

NV 

I 

Number  of  data  points  per  record,  provided  by 
an  information  sheet  called  DIGITAL  TAPE  LOG. 

CWA 

R 

Cell  constant,  set  to  0.505?. 

CW1 

R 

Cell  width  in  microsecs,  provided  by  DIGITAL 

TAPE  LOG. 

NA 

I 

Number  of  cells  integrated,  provided  by  DIGITAL 

TAPE  LOG. 

NORM 

I 

Range  normalization  control,  set  to  0.0  for  no 
range  normalized  data,  set  to  nonzero  value  for 
range  normalized  data,  provided  by  DIGITAL  TAPE 

LOG . 

MINV 

I 

Minimum  data  value  (bit  count)  for  noise  level 
cutoff  from  calibration  curve. 

MAXV 

I 

Maximum  data  value  (bit  count)  above  which  all 
data  is  in  error  from  calibration  curve,  normally 
set  to  255. 

I CUT OF 

I 

Data  value  (bit  count)  below  which  leftmost  portion 
of  calibration  curve  begins  to  deviate  from  remainder 
of  curve. 

FA 

R 

Intercept  of  calibration  curve  on  DBM  axis. 

GLC 


LMBDA 

AB 

BB 

RAINHT 


FALEFT 


R 

FM  R 

FMLEFT  R 

NVAL  R 


Intercept  of  leftmost  part  of  calibration  curve  on 
DBM  axis,  add  +115  DBM  to  plot  DBM  values  to  obtain 
correct  intercept  value. 

Slope  of  calibration  curve. 

Slope  of  leftmost  part  of  calibration  curve,  slope 
defined  as  :  DELTA  DBM/DELTA  BIT  COUNT. 

Currently  not  used,  set  to  60*0.0. 


3.  TIME  AND  MEMORY  REQUIREMENT 


The  1 07K 
Approx imately 
file. 


memory  in  octal  is 
40  cp  seconds  processing 


requi red 
time  is 


to  run  this  program, 
required  for  each 


4.  TAPER 

TAPE1  is  an  input  tape,  it  is  a  binary  tape,  odd  parity.  There 
is  more  than  one  file  per  tape,  each  file  has  to  be  processed 
separately.  The  deck  is  set  up  to  copy  TAPE1  to  the  tape  to  skip 
bad  files.  TAPE 5  is  for  punched  input  data.  TAPER  is  for  printed 
output,  the  deck  is  set  up  to  rewind  TAPER  and  copy  it  to  the 
output  file.  The  structure  and  information  of  the  input  tape  is 
provided  by  a  sheet  called  the  DIGITAL  TAPE  LOG.  TABLE  A-3  is  a 
sample  sheet  of  the  DIGITAL  TAPE  LOG. 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REOUIRED 

ATNCOR  (NEWVLC ,CL, S INEL , RNG , DY , YY ) 

CONVRT  (IVAL, NEWVLC) 

READW  ( I  EOF) 

RNCGEN  ( CW 1 , NA , I SAVE ) 
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PROGRAM  IWCCI 


1.  DESCRIPTION 

Program  IWCCI  was  modified  from  program  LWC  which  was  used  in 
the  AFGL-2  USSR  eleven  station  analysis.  Both  programs  were 
developed  to  read  AFGL-2  meteorological  data  cards.  Program  IWCCI 
classifies  Cirrus  clouds  according  to  their  amount,  from  0  to  10 
tenths  coverage;  by  type,  whether  single  or  multi-layer  and  whether 
associated  with  Cumulonimbus  clouds  or  not;  and  also  lists  their 
base,  top,  liquid  water  content  and  thickness  (top-base) .  It  also 
produces  a  statistical  table  in  which  the  monthly  number  of 
occurrences  by  time  (each  three  hour  period  of  day)  and  by  sky 
cover  amount  (0  to  10  tenths)  and  the  summation  of  each  is  listed. 
The  total  number  of  clear  cases  is  also  tabulated  by  month  and 
listed  beneath  the  table. 

Input  to  the  IWCCI  program  consists  of  the  meteorological  data 
cards  generated  from  the  AFGL-2  eleven  station  analysis  or  similar 
analysis. 

Output  from  the  IWCCI  program  consists  of  a  printout  comprised 
of  twelve  tables.  Tables  one  to  eleven  give  Cirrus  cloud 
information  based  on  Cirrus  cloud  sky  coverage  from  0  to  10  tenths 
and  the  twelfth  table  lists  the  total  monthly  occurrance  of  each 
Cirrus  event. 

2.  DECK  SET  UP 

7.1  CONTROL  CARDS  (IWCCI  in  UPDATE  FILE  under  TUNG  ID) 

JQBNM ,CMJ05000,T20.  ACT  NO.  NAME 
ATTACH , OLDPL , TWCC I , TD=TUNG . 

UPDATE, N. 

FTN ,  I=C0M.PILE,  SL,  PL  =  77 77777. 

LDSET (PRESET=ZERO) 

LGO. 

7/8/9 

* IDENT  CHECK 
♦DELETE  IWCCI. 29 
MTH=1 
*C  IWCCI 
7/8/9 

METEOROLOGICAL  DATA  CARDS 

f/vr./o 

(MTH  represents  the  month  of  data  being  processed  with  1=JAN, 
?=FFP,  etc..  Only  one  month  of  data  may  be  processed  at  a  time  and 
o n  1  v  one  s':,-*-  ion  ray  be  processed  at  a  time.  The  month  designated 
by  MTH  rust  be  ‘he  sane  as  the  month  of  met  eo  r  o  1  og  i  ca  1  data  cards 
be  i  no  ]  r ore ssr d . ) 


/  •  f 


’NPUT  data  cards 


The  input  data  cards  consist  of  the  meteorological  data  cards. 
The  format  for  the  meteorological  data  cards  is  given  in  TABLE  A-4. 

3.  TIME  AND  MEMORY  REQUIREMENT 

The  105k  in  octal  of  memory  is  required  to  run  this  program. 
Approximately  10  cp  seconds  processing  time  is  required  to  process 
one  month  of  data  for  one  station. 

4.  TAPES 

TAPE3  is  used  to  store  ESI  profiles  which  is  not  used  in  this 
version  of  the  program.  TAPE5  is  reserved  for  the  input  data  cards 
(format  in  TABLE  A-4).  TAPE 6  to  TAPE8  are  reserved  for  the 
statistics  tables  output  which  are  not  used  in  this  version.  TAPE11 
to  TAPE 21  are  used  to  store  and  output  eleven  sets  of  data 
classified  by  sky  coverage. 

5.  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

REARAY (Y1 , Y?,XX, Y4, Y5,Y*,KK) 

ZER0(Y1,Y2, XX, Y4,Y5, Y6.KK) 

STATTS (NO, IJK) 

MEDIAN ( I , J) 

INT7.  (B,L,X,Y,21  ,7.2,7.) 
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TABLE  A- 1 


AFGL  ALCOR  WEATHER  TAPE 

The  format  and  contents  of  the  AFGL  ALCOR  Amplitude  or  Phase 
weather  tape  to  be  generated  on  the  CDC-6600  is  described  below. 
The  tape  will  be  7-track,  BOO  BPI,  Fortran  unformatted. 


The  first  record  on  the  tape  will  be  a  header  record  and  all 
others  will  be  data  records.  All  values  are  store  in  integer  form, 
a  factor  is  provided  for  division. 

HEADER  RECORD 


NAME 

TYPE 

MCON  ARRAY 

WORD 

DESCRIPTION 

IDESCR ( 1 ) 

EBCDIC 

1 

1 

40  characters  from  input  card 

IDESCR (10) 

10 

10 

ITBAND 

FX 

11 

1 1 

=0  (NB) ,  =1  (WB) 

ITREEL 

FX 

12 

12 

Reel  number 

IMTH 

FX 

1  3 

13 

Month 

IDAY 

FX 

14 

14 

Day 

IYR 

FX 

15 

15 

Year  (last  two  digits) 

ILNCH 

FX 

16 

16 

Lift  off  (mi 1 1 i secs) 

IPULST 

FX 

17 

17 

First  pulse  No.  on  tape(INPUT) 

IPULSP 

FX 

18 

18 

Last  pulse  NO.  on  tape(INPUT) 

TRKGTE 

FX 

19 

19 

ALCOR  tracking  gate  *100 

GTESPC 

FX 

20 

20 

Tntergate  spacing  (m)  *100 

- 

21- 

21- 

Spares 

DATA 

RECORD 

NAME 

TYPE 

NCON  ARRAY 

WORD 

DESCRIPTION 

RCSLC  (1  ) 

FX 

1 

1 

LC  Amplitude  or  Phase  (db)  *100 

RCSLC (?) 

FX 

2 

1 

For  Amplitude, if  no  sign  assume 

RCSLC ( 3 ) 

FX 

•? 

1 

negtive  numbers. 

RCSLC (4) 

FX 

4 

1 

For  Phase, the  left-most  bit  of 

— 

- 

- 

the  15  bit  integer  is  sign  bit. 

RCSLC (169) 

FX 

If  9 

43 

RCSLC (170) 

FX 

170 

43 

RCSRC(l) 

FX 

171 

43 

RC  Amplitude  or  Phase  (db)  *100 

RCSRC ( 2 ) 

FX 

172 

4  3 

RCSRC (170) 

FX 

340 

85 

GMT 

FX 

341 

86 

GMT  (Seconds)  *1000 

_ 

342 

86 

TAL 

FX 

343 

86 

TAL  (seconds)  *1000 

- 

344 

86 

NUMPRI 

FX 

345 

87 

Pul se  number 

IHRS 

FX 

346 

87 

Hours 
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NAME  TYPE  NCON  ARRAY  WORD  DESCRIPTION 


IM  IN 

FX 

347 

87 

I  SEC 

FX 

348 

87 

IFRAC 

FX 

349 

88 

RANGE 

FX 

350 

88 

AZ 

FX 

351 

88 

EL 

FX 

352 

88 

CALT 

FX 

353 

89 

- 

FX 

- 

- 

- 

FX 

400 

100 

Minutes 
Seconds 
Mill  iseconds 

range  to  tracking  gate  (m)  *10 

Azimuth  (deg)  *100 

Elevation  (deg)  *100 

Altitude  (m)  *10 

Spares 

Spare 


NOTE  :  In  data  records,  the  array  NCON  are  packed  four  to  each  CDC  word. 
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TABLE  A-? 


AFGL  SPECIAL  TRADEX  WEATHER  TAPE 


cal  ibrat  c°ntents  of  a  special  TRADEX  S-narrow  band  chip 

b:i^e.tatp:Peto  ^igenKati^a^r0r°  a/g^ 

unformatted,  containing  standard  SCOPE  records.  '  °r  ran 


of.h  JirKfc  5eC°rd  °n  the  tape  wiU  be  3  header  record  and  all 

others  will  be  data  records.  The  header  record  will  contain 

information  describing  the  contents  of  the  tape.  The  data  record 

:  1  a  ^  ^retriC.lnf°rati0n  updated  every  s  and  one  pulse 

or  calibrated  signature  data. 


HEADER  RECORD  (20  WORDS) 


WORD 

1 

2 

3 

4-10 

11 

1? 

1  3 

14-20 


DESCRIPTION 


Polarization  (9=PP,  10=OP) 
No.  of  range  gates  per  pulse 
Intergate  spacing  (m) 

Tape  description 

Tape  start  (GMT  total  s) 

Tape  stop  (GMT  total  s) 

Gate  number  of  tracking  gate 
Spares 


TYPE 

Integer 

Integer 

Real 

Display  code 
Integer, LSB  is  0.0001 
Integer, LSB  is  0.0001 
In  teger 


DATA  RECORD 


WORD 


6-1  0 
11- 


DESCRIPTION 


Pulse  time  (GMT  total  s) 

Range  (m) 

Azimuth  (deg) 

Elevation  (deg) 

Time  (GMT  total  s)  associated 

with  R, AZ ,and  EL 

Spares 

Rader  cross  section  for  all 
gates  of  the  pulse  packed  four 
to  each  CDC-K600  word. 

The  left-most  bit  of  the  15  bit 
integer  is  the  sign  bit. 
Nagative  integers  are  in  one's 
complement  form. 


TYPE 

Integer, LBS  is  0.0001 

Real 

Real 

Real 

Integer, LBS  IS  0-0001 


Integer, LSB  is  0.1  dbsm 
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TABLE  A— 4 


COLUMNS 

1 

2-4 

5-r, 

•’_o 
in-]  i 

]  ?-n 

14-17 

ip-21 

22 

23 


24-27 

28-3  3 

32-33 

34-35 


3  r,  -  a.  p 

4  1  -58 
59-7S 

NOTE: 

If 

prof  il  e, 
all  data 


vFTEOROLOG  rCAL  DATA  CARD  FORMAT 
SAMPLE  DATA 


]  55 

0  1 
FEB 
73 
00 

0200 
03  20 


2 


01.4 
05.0 
1  0 

type 


0.100 


Blank. 

Location  (station  number, 

15  5  =  Mo  scow 1  . 

Day. 

Month. 

Year. 

Time  (007,037, etc.). 

Area  representativeness 

(200  nm  in  the  N-S  direction). 

Area  representativeness 

(120  nm  in  the  F-W  direction). 

Blank  . 

Identifier  (2  =  precip.  water 
content  in  g/m**?;  ]  =  liquid  or 
cloud  water  content  in  g/m**3; 
blank  =  clear)  . 

Altitude  (lower)  of  layer  in 
thousands  of  feet  (1400  feet). 
Altitude  (upper)  of  layer  in 
thousands  of  feet  (5000  feet). 
Coverage  in  tenths  (10/10). 

Type  of  precipitation  or  cloud 
(T.C.  =  ice  crystal; 

S.S.  =  small  snow; 

L.S.  =  large  snow; 

RN.  =  rain  with 

cloud  type  symbols  those 

normal  1 y  used)  . 

Water  content  in  g/m**1. 

Second  and  third  sets  of  data; 
follow  same  format  as  23-40. 


more  than  three  sets  of  data  are  present  during  one 
columns  1-22  are  repeated  on  the  following  card(s)  until 
i s  tabu) ated . 
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APPENDIX  B 

THE  WEATHER  CORRELATION  RUNS  OF 
26  JUNE  1979  AT  KWAJALEIN 
David  W.  Blood 
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APPENDIX  B 


THE  WEATHER  CORRELATION  RUNS  OF 
26  JUNE  1979  AT  KWAJALEIN 


1.  Introduction 

The  weather  correlation  test  series,  KMR  4028,  was  conducted  at 
Kwajalein  Atoll,  Kwajalein  Missile  Range  on  26  June  1979.  The  purpose 
of  the  linked- radar/aircraft  experiment  was  to  collect  data  simultaneously, 
both  by  radar  backscatter  and  by  on-bcard  aircraft  measurements  of  the 
hydrometeor  type  and  size  distribution  for  correlation  and  calibration 
purposes.  The  tests  were  conducted  to  establish  liquid  water  versus 
reflectivity  relationships  for  a  re-entry  mission  where  radar-only  data 
would  be  available.  The  series  consisted  of  several  runs  (numbered  1 
through  5)  during  the  period  of  0200  and  0400  C.MT  on  the  GMT  date  indicated 
where  the  Lear-jet  aircraft  flew  different  passes  through  weather  storm 
cells  while  being  tracked  by  an  MPS-36  radar  on  Kwajalein,  Is.  at  the  Atoll. 
The  Lincoln  KREMS  radars  (ALCOR,  C-band  and  TRADEX,  S-band)  were  linked  to 
the  tracking  radar  and  designated  to  collect  backscatter  returns  from  a 
cell  volume  immediately  (3  km)  ahead  of  the  aircraft.  The  radar  data  were 
later  processed  off-line  by  a  program  called  "moist"  to  calculate  quantities 
related  to  atmospheric  moisture  content.  This  data,  in  turn,  may  be  cor¬ 
related  against  the  aircraft  on-board  measurements  to  check  the  calibra¬ 
tion  constants  and  procedures  used  in  converting  radar  backscattering 
cross  sections  to  radar  reflectivity  factor,  Z,  which  in  turn  is  relatable 
to  the  moisture  content  within  the  range  cell  of  question.  Ideally,  the 
two  backscattering  radars,  having  roughly  the  same  resolution  (with  the 
waveforms  used  for  this  type  of  measurement),  will  have  agreement  to 
within  about  I  dB  of  each  other  in  terms  of  the  backscattered  Z  levels. 

On  this  series,  however,  considerable  disagreement  was  experienced  on 
some  of  the  runs  and  in  certain  regions  of  a  particular  run  amounting  to 
differences  at  times  up  to  25  dB  between  the  two  radars.  This  rendered 
questionable  the  use  of  either  radar  measurements  at  these  times.  In 
this  report,  only  the  radar  operations  and  results  will  be  discussed. 

The  operations  consisted  of  four  runs  in  the  vicinity  south  and 
slightly  east  of  Kwajalein,  Is.  at  altitudes  of  4500  to  4800  meters  and 
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one  run  north-east  of  Roi-Namur,  Is.  at  altitudes  of  2200  to  2600  meters. 
Correspondingly,  the  KREMS  radars  were  directed  to  look  toward  the  south¬ 
east  at  ranges  varying  from  75  to  115  kms  (elevations  from  2.9  to  1.8 
degrees)  for  run  nos.  1  to  4  and  toward  the  northeast  at  ranges  varying 
from  15  to  22  kms  (elevations  from  8.5  to  6  degrees)  on  run  no.  5. 

Figure  1  shows  the  location  of  the  tracks  in  terms  of  radar  coordinates 
(slant  range  and  azimuth)  from  Roi-Namur.  It  appears  that  for  the  most 
part,  the  radar  operations  were  nominal  in  the  linked  mode  in  following 
the  designated  tracks  and  in  monitoring  data  within  the  dynamic  range  of 
the  radar  receivers  through  the  real-time  attenuation  switching  algorithms. 
There  were  certain  losses  of  tracking  (and  of  maintaining  locked  radars) 
at  the  mid-portion  of  run  no.  1,  for  example,  and  at  the  beginning  and 
end  of  run  no.  5  (dashed  traces  on  Figure  1).  Since  these  regions  may 
easily  be  edited  out  of  the  acceptable  data  sets  no  further  discussion 
will  be  devoted  to  these  limited  difficulties. 

The  TRADF.X  S-band  radar  antenna  system  was  operating  in  a  45°  linear 
polarization  configuration  as  opposed  to  the  circular  polarization  con¬ 
figuration  previously  used  and  also  used  on  ALCOR.  This  caused  no  unusual 
difficulties  in  terms  of  calibration  or  antenna  patterns,  as  initially 
suspected  for  the  radar  disagreements.  The  goodness  of  radar  calibrations 
was  verified  by  (20  inch)  sphere  calibrations  conducted  one  day  after  the 
weather  correlation  runs.  Both  TRADF.X  and  ALCOR  obtained  sphere  back- 
scattered  RCS  measurements  agreeing  with  the  theoretical  optical  RCS 
value  of  -7.0  dBsm  to  within  ±1  dB.  The  polarization  isolation  between 
the  principal  (PP)  and  opposite  (OP)  polarization  received  channels  on 
these  high  signal-to-noise  (S/N)  checks  exceeded  27  dB  on  both  radars. 
Antenna  pattern  checks  with  the  linearly  polarized  TRADEX  S-band  config¬ 
uration  were  requested  and  they  indicated  no  change  from  the  earlier  cir¬ 
cular  configuration  patterns.  The  question  remained,  therefore,  as  to 
why  the  large  discrepancy  in  backscattered  dBZ  values  between  the  two 
radars  supposedly  trained  upon  the  same  weather  cell. 


Figure  1  Weather  correlation  runs,  June  26,  1979 


2.  Observations  on  June  26,  1979 


Figures  2  through  6  show  the  radar  measurements  in  dBZ  for  runs  no.  1 
through  5  respectively.  The  ALCOR  measurements  (dashed)  may  be  compared 
with  TRADEX  (solid)  for  each  of  these  runs.  It  is  evident  that  a  large 
departure  occurred  between  TRADEX  and  ALCOR  during  all  portions  of  runs 
1  and  2  and  during  the  latter  portion  only  on  run  3.  Corresponding  to 
the  dBZ  plots,  the  received  polarization  ratio  (principal  to  opposite, 

PP-OP  in  dB)  is  shown  at  the  bottom  of  each  figure.  During  the  periods 
of  backscatter  disagreement  the  polarization  ratios  were  also  generally 
poorer  than  would  be  expected  for  rain  scatter. 

As  a  result  of  the  observed  polarization  ratio  degradation  and  an 
examination  of  the  theoretical  ratio  to  be  expected  from  rain, a  set  of 
editing  rules  were  developed  for  selecting  out  the  usable  data  as  indicated 

in  Table  1 . 

Several  other  observations  may  be  made  from  the  radar  data.  The  early 
portion  of  run  no.  3,  run  no.  4  and  run  no.  5  all  show  favorable  back¬ 
scatter  agreement  (to  within  6  dB)  between  TRADEX  and  ALCOR.  At  these 
times  the  polarization  ratios  were  generally  high  (exceeding  18  dB)  when 
signal-to-noise  conditions  were  high.  When  agreement  was  favorable  on 
runs  3  and  4,  there  still  remained  a  6  dB  difference  between  TRADEX  and 
ALCOR,  ALCOR  being  lower.  Finally,  on  run  5,  much  nearer  the  KREMS  radars 
the  dBZ  levels  agree  very  well  to  within  ±  a  few  dB  as  has  been  observed 
on  other  normal  missions. 

3.  Possible  Explanation  of  Results 

It  is  believed  that  the  observed  different  dBZ  levels  and  PP/OP 
degradation  are  not  the  result  of  a  single  cause  but  a  combination  of 
effects.  The  effects  all  may  be  classified  as  propagation  effects,  however, 
as  opposed  to  hardware  or  tracking  problems.  There  are  several  reasons 
to  believe  that  the  radars  were  operating  correctly.  The  sphere  calibra¬ 
tions  cited  earlier  strongly  suggest  no  major  biases  in  the  radar  set-up. 
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•  M402I  CO**  *U*  | 


(— )TRADEX 
( — )  ALCOR 


HB)  DBZ  VALUES 

34.0  -30*0  -14.0  2.0  18.0  34.0  SO-O 


PP-OP  (dB)  OBZ  VRLUES 

?.n  18.0  34.0  -30.0  -14.0  2.0  18*0  34.0  50. 


TABLE  1 


MINIMUM  EDITING  RULES  FOR  SELECTING 
USABLE  RADAR  BACKSCATTER  DATA  FROM  RAIN 


Range  of  dBZ 


Minimum  Rule* 


Ideal  (Theory) 


If  dBZ  <  0, 

then 

PP-OP 

> 

20 

dB 

off  S' 

0  £  dBZ  <  10, 

then 

PP-OP 

> 

19 

dB 

- 

10  £  dBZ  <  20, 

then 

PP-OP 

> 

18 

dB 

30 

20  £  dBZ  <  30, 

then 

PP-OP 

> 

17 

dB 

25 

30  £  dBZ  <  40, 

then 

PP-OP 

> 

15 

dB 

20 

40  £  dBZ  <  SO, 

then 

PP-OP 

> 

13 

dB 

17 

SO  £  dBZ  <  60, 

then 

PP-OP 

> 

12 

dB 

15 

60  £  dBZ 

then 

PP-OP 

> 

11 

dB 

14 

♦Somewhat  arbitrary  at  lower  dBZ  values  due  to  lack  of  knowledge  of 
noise  statistics  influencing  the  PP /OP  ratio. 
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Some  runs  were  good  and  some  went  from  good  toward  a  worsening  condition 
(run  3  for  example).  In  no  case  when  radar  attenuators  were  switch-in, 
thus  maintaining  the  proper  dynamic  range,  did  "step  function"  changes 
occur  in  the  dBZ  values;  this  suggests  that  there  was  never  an  improper 
calibration  constant  in  use. 

On  run  3  when  conditions  worsened,  the  region  affected  was  localized. 
Figure  1  shows  darkened  regions  along  the  trajectory  of  the  tracked 
scattering  volume  which  correspond  to  high  regions  of  backscatter.  The 
beginning  and  end  times  of  these  scatter  regions  permit  a  crude  mapping 
of  the  boundaries  of  the  storm  regions  outlined  on  the  figure.  The  region 
of  poor  backscatter  agreement  is  also  shaded  in  Figure  1.  Upon  closer 
examination  of  the  azimuthal  weather  scans  taken  at  ALCOR  at  0135  and 
0334  GMT  at  low  elevation  angles  (2.3°)  a  "clutter  notch"  was  observed  to 
exist  at  precisely  the  azimuths  in  question  (140°  to  155°  azimuth).  This 
notch  can  account  for  at  least  10  dB  in  the  drop  of  ALCOR  dBZ  levels  when 
looking  at  these  azimuths  at  low  elevation  angles.  It  should  be  noted 
that  the  radar  elevation  angles  being  designated  to  the  KREMS  radars 
ranges  from  2.0  to  2.5  degrees  in  the  shaded  region  of  Figure  1.  Figure  7 
is  a  tracing  of  the  azimuthal  scan  taken  at  ALCOR  at  the  conclusion  of 
run  no.  4  (2.3°  elevation).  A  "clutter  notch"  is  shown  at  the  indicated 
azimuths.  Correspondingly,  the  rain  scatter  region  of  interest  (R  =  80  to 
110  km]  shows  a  reduced  amplitude  (speckled  pattern)  versus  the  more  solid 
region  (cross  hatched)  observed  at  azimuths  greater  than  155°.  It  is  this 
localized  "attenuation"  effect  which  is  the  prime  cause  of  the  lower 
ALCOR  scatter  observed  in  Figures  2,  3  and  4.  The  cause  of  this  is 
probably  attributable  to  local  palm  tree  obscuration  (screening)  at  these 
particular  azimuths  at  low  elevations  affecting  only  the  ALCOR  radar, 
physically  set  into  and  behind  the  trees. 

The  local  screening  effect  at  ALCOR  does  not  adequately  explain 
the  poor  PP/OP  ratios  at  both  TRADEX  and  ALCOR  during  these  times. 

Another  feature  identified  in  the  0334  GMT  scan  at  ALCOR  (Figure  7)  was 
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Figure  7  ALCOR  azimuthal  weather  scan  -  2.3°  elevation 


that  classed  as  a  "thin  layer  return"  over  similar  azimuthal  regions. 

On  the  corresponding  vertical  scan  taken  at  0339  (not  shown,  179°  azimuth) 
this  return  may  be  identified  as  a  thin  (in  vertical  extent)  layer  type 
of  return  perhaps  related  to  a  refractive  duct  being  present  producing  an 
apparent  2  km  high  return  mid-way  between  the  radars  and  the  scattering 
volume  of  interest.  It  is  not  known  whether  this  is  an  actual  backscattered 
return  at  this  height  or  perhaps  more  likely  a  "ghost"  return  of  sea 
clutter  produced  at  the  proper  elevation  angle  by  rays  being  "refracted- in" 
by  a  ducting  layer.  Indeed,  evidence  was  available  from  an  earlier  radio¬ 
sonde  profile  that  a  strong  ducting  condition  existed  in  the  refractivity 
versus  height  profile  (Figure  8).  If  the  ducting  layer  prevailed  at  the 
indicated  height  over  the  Atoll  it  would  have  refracted  rays  at  ranges  of 
only  15  to  25  km  producing  the  apparent  sea-scatter  return  at  the  30  to  50 
km.  In  any  case,  the  thin  layer  return  of  Figure  7  may  be  an  indicator 
of  refractive  conditions  which  could  be  responsible  for  the  polarization 
ratio  degradation  on  TRADEX  and  ALCOR  by  producing  a  tropospheric  multi- 
path  type  of  scattering  return  from  the  rain  volume  of  interest.  The 
latter  explanation,  though  somewhat  tenuous  due  to  inadequate  supporting 
data,  does  point  to  the  hazards  of  low  elevation  measurements  and  of 
collecting  data  during  periods  of  potential  duct  situations. 

An  analysis  of  the  ALCOR  extended  range  gate  data  (about  the  tracking 
cell  range)  was  performed  in  search  of  possible  time  sidelobe  interference 
produced  by  a  hard  return  from  the  aircraft  itself.  All  of  the  poor  dBZ 
agreement  regions  were  processed  from  the  ALCOR  narrow  band  waveform  data. 
The  result  from  this  study  was  that  no  evidence  of  a  hard  body  (discrete 
return)  target  was  identified  as  affecting  the  ALCOR  radar  data. 

4.  Conclusions 

1.  It  appears  that  the  ALCOR  data  taken  on  runs  1,  2  and  parts  of  3 
are  unusually  low  in  amplitude  due  to  a  local  screening  problem  identified 
in  the  radar  sea  clutter  return. 
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Figure  8  Refractivity  profile  at  Kwaialein  for  2300  GMT  on  25  June  1979 
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2.  The  lower  PP/OP  ratio  on  these  runs  appears  to  be  due  to  a 
totally  separate  propagation  problem  affecting  both  TRADEX  and  ALCOR 
at  the  low  elevation  angles  0  <  2.5°  and  long  ranges  R  >  80  km.  This 
may  be  attributable  to  a  refractive  layer  effect  over  the  Atoll. 

3.  The  minimum  editing  rules  developed  (Table  1)  should  permit  a 
technique  for  selecting  useful  rain  backscatter  data  for  weather  correla¬ 
tion  analysis. 

5.  Recommendations  for  Next  Series 

1.  Low  elevation  angle  ALCOR  data  should  be  avoided  due  to  local 
terrain  features  at  certain  azimuths  from  Roi-Namur. 

2.  Low  altitude,  long  slant  range  runs  should  be  avoided  when 
possible  with  KREMS  or  the  tracking  radar  unless  radiosonde  measurements 
show  no  evidence  of  refractive  layers  present  at  the  time  of  the  weather 
runs. 

3.  Sphere  calibrations  are  advisable  in  the  immediate  time  frame 
of  the  tests. 

4.  In  general,  if  only  one  radar  were  available,  TRADEX  would  be 
preferred  over  ALCOR. 
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APPENDIX  C 


CIRRUS  CLOUD  PARTICLE  DETECTABiLITY 
USING  THE  KWAJALEIN  RADARS 


1.  Introduction 

Recent  investigations  of  clear  air  turbulence  backscatter  with  the 
Kwaialein  (KREMS)  radars  has  led  to  a  reassessment  of  cirrus  cloud  (ice 
particle)  detectability  by  making  use  of  coherent  signal  processing  tech¬ 
niques.  The  ability  to  make  measurements  of  low  particle  densities  of 
tenuous  clouds  by  radar  is  of  interest  to  the  question  of  re-entry 
vehicle  (RY)  surface  erosion  and  the  triggering  of  the  transition  between 
laminar  and  turbulent  flow  which  may  occur  at  altitudes  of  about  10  km 
for  otherwise  smooth,  slender  re-entry  bodies.  The  particle  sizes  of 
concern  are  ice  crystals  of  diameter  greater  than  70  microns  and  at  densi¬ 
ties  of  only  a  few  particles  per  cubic  meter.  At  these  levels  of  occur¬ 
rence,  sensitive  radars  may  provide  the  practical  means  to  map  remotely 
the  atmospheric  environment  along  the  re-entry  corridor. 

For  clear  air  turbulence  (CAT)  backscatter  measurements,  the  radial 

component  of  the  wind  may  be  sensed  due  to  the  Doppler  shift  of  the  CAT 

scattered  energy^.  The  same  mechanism  may  be  anticipated  for  cloud 

particles  in  that  the  spectral  returns  for  radar  viewing  angles  off- 

vertical  will  be  Doppler  shifted  by  the  prevailing  winds  at  altitude. 

If  viewing  azimuths  are  other  than  tangential  to  the  wind  vector  (i.e. 

toward  or  away  from  the  wind  direction)  the  received  signals  are  shifted 

spectrally  away  from  the  stronger  sea  and  land  clutter  signals  obtained 

near  zero  Hertz  (0  m/s  velocity).  The  coherent  processing  capability 

therefore  affords  the  possibility  of  several  tens  of  dB  improvement  over 

2 

the  non-coherent  detectable  level  of  the  radars  previously  considered". 

The  required  geometry  is  achieved  at  the  Kwaialein  Atoll  with  KREMS 
radars  viewing  the  RV's  toward  the  prevailing  winds  for  typical  situa¬ 
tions.  Figure  l  shows  an  example  of  CAT  returns  sensed  using  the  TRADEX 
(I.-band)  radar  at  a  range  corresponding  to  the  height  of  5.7  km  (Ref.  3). 
A  similar  return  to  the  "turbulence"  return  may  be  expected  from  ice  par¬ 
ticles  at  higher  altitudes  and  with  radar  velocity  components  correspond¬ 
ing  to  the  winds  at  cloud  heights. 
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F  igure  I  Sample  velocity  spectrum  obtained  from  TRAD1.X 


C-3 


2.  Weak  Cirrus/Ice  Detectability 


With  the  ALTAIR  (VHF/UHF)  and  TRADEX  (L  5  S  band)  radars  several 
high-energy  waveforms  exist  which  have  Doppler  processing  potential 
together  with  adequate  range  resolution  for  consideration  in  cirrus  cloud 
particle  detection.  The  unambiguous  velocity  window  (for  the  radar  wave¬ 
lengths)  ranges  from  40  to  60  m/s  along  the  radar  line-of-sight.  This 
window  has  proven  adequate  for  most  wind-measurement  observations.  The 
ALCOR  (C  band)  radar,  although  having  a  smaller  wavelength  appropriate 
for  ice  crystal  detection,  has  limitations  in  energy  and  unambiguous 
velocity  (a>  6  m/s)  which  do  not  provide  an  overwhelming  advantage  over 
the  use  of  the  lower  frequencies.  Table  1  lists  the  radar  parameters 
and  comparative  detectabilities  of  the  more  appropriate  waveforms.  From 
the  list,  the  TRADF.X  L-chirp,  S-chirp  or  ALTAIR  UCWL  waveforms,  for 
example,  should  provide  alternatives  to  the  use  of  C-band  and  make  the 
separation  in  the  spectral  domain  more  feasible. 

2 

The  minimum  detectable  structure  constant  Cn*-  i?  computed  in  Table  1 
at  a  5  db  S/N  level  and  for  a  radar  slant  range  of  42  km.  This  range 
corresponds  to  a  height  of  about  10  km  but  within  the  re-entry  corridor 
for  low  angle  target  trajectories  impacting  short  of  the  lagoon  from 
Vandenberg  AFB.  The  length  of  observational  time  used  throughout  is  80 
seconds,  the  duration  found  adequate  for  CAT  backscatter/wind  measure¬ 
ments.  Equivalently,  the  reflectivity  factor  Z  and  the  structure  con- 
2 

stant  Cn  are  related  by: 


Z 


7  1 R 

0.3786  (Cn  x  10  )  X 

Ik*  I2 


11/3 


(mm^/m3) 


(1) 


where  X  =  radar  wavelength  (m) , 

2  -2/3 

Cn  =  structure  constant  (m  '),  and 


=  0.209  for  ice  particles. 


It  is  evident  that  5  dB  S/N  levels  can  be  obtained  for  dBZ  values  as  low 
as  -50  dBZ  with  coherent  processing  and  at  the  42  km  slant  range.  It  is 
necessary  for  the  radar  beam  to  dwell  for  an  80-second  observational 
period  at  the  desired  azimuth  and  elevation  angle  to  obtain  these  sensi¬ 
tivities.  Since  it  is  only  the  detection  of  particles,  not  the  Doppler 
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velocity  that  is  of  interest,  it  is  not  necessary  fov  the  radar  to  scan 
in  azimuth  as  was  the  requirement  for  wind  measurements. 

The  ice  particle  densities  are  also  listed  in  Table  1  for  ice 
particles  greater  than  70  and  80  microns  in  diameter  as  derived  from 
Burke  et  al.,  1978.  Possible  waveforms  for  detecting  tenuous  particles 
(lower  densities  of  a  few  per  meter^l  would  be  the  ALTAIR  (UHF)  UCWL 
waveform  or  TRAOEX  L  or  S  band  chirp  waveforms.  The  ALTAIR  UXL  and 
VCWl,  waveforms,  though  nigh  energy,  may  not  offer  adequate  height  resolu¬ 
tion  (1.45  km  and  1.09  km  respectively  at  the  radar  range  resolutions  of 
6000  meters  at  UXL  and  4500  meters  at  VCWL)  for  the  measurements  at  42  km 
slant  range. 

With  the  TRADRX  L-chirp  waveform,  the  lowest  particle  densities  are 
obtained.  Forty-two  particles  per  cubic  meter  should  be  detectable  for 
particle  diameters  larger  than  80  microns.  Correspondingly,  the  same 
waveform  could  detect  only  three  particles/m'^  at  near  overhead  geometries 

3 

with  a  5  dB  S/N  or  15  particles/m  at  42  km  range  but  with  unity  S/N 
ratio.  A  longei  integration  time  could  further  improve  upon  these  sensi¬ 
tivities.  The  corresponding  altitude  resolution  would  be  near  200  meters 
at  high  elevation  angles  or  50  meters  at  42  km  range. 
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3.  Operational  Configuration 


It  will  be  necessary  to  distinguish  between  clear  air  turbulence 
returns  (in  the  absence  of  ice)  and  scatter  from  cirrus  cloud  ice  par¬ 
ticles  alone.  This  may  be  accomplished  by  the  use  of  two  frequencies 

illuminating  the  same  volume.  If  CAT  provides  the  scatter  mechanism, 

2 

the  two  radars  should  sense  the  same  structure  parameter,  Cn  at  the 

same  altitude.  Figure  2  shows  altitudes  where  cirrus  particle  back- 

scatter  could  not  be  easily  distinguished  on  a  single  radar  frequency 

from  CAT  returns  (shown  from  Ref.  3)  because  of  similar  levels  of  back- 

scattered  returns.  Given  only  cirrus  particles  at  a  particular  altitude, 

2 

a  different  equivalent  Cn  level  would  be  sensed  with  two  frequencies 

and  the  more  sensitive  radar  would  dominate  in  detecting  these  particles. 

It  should  be  possible  therefore  to  eliminate  CAT  as  the  scatter  mechanism 
2 

in  question  on  a  Cn  display  containing  two  radar  frequency  outputs 
versus  altitude.  A  three -frequency  (two  radar)  technique  would  provide 
even  a  more  positive  identification  of  scatter  type. 


4. 


Conclusions 


The  sensitivity  calculations  of  Table  1  indicate  the  potential  of 
KREMS  radars  in  detecting  cirrus  cloud  ice  particles  at  low  levels  of 
density  within  the  re-entry  corridor.  Alternatives  for  even  more  sensi¬ 
tivity  are  to  search  at  higher  elevation  angles  or  integrate  for  periods 
longer  than  80  seconds.  Coherent  (Doppler)  processing  of  the  returns  to 
gain  this  sensitivity  is  most  promising  by  using  the  TRADEX  and  ALTAIR 
waveforms  as  opposed  to  the  ALCOR  C-NB  waveform.  A  two  (or  three) 
frequency  method  should  offer  a  means  of  positive  identification  of  the 
particle  returns  from  those  produced  by  turbulence.  Using  the  multiple 
frequency  observation  technique,  80  urn  particles  with  densities  as  low 
as  42  m-3  would  be  detectable  (TRADEX  L-chirp)  at  a  range  of  42  km  in 
the  presence  of  turbulence  of  the  intensities  observed  to  date  at 
Kwajalein. 
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APPENDIX  D 


THE  AFGL-1  DECISION  TREE  FOR 
LIQUID  WATER  CONTENT  ASSIGNMENT 


1.  Introduction 

The  AFGL-1  model  is  a  technique  for  converting  standard  meteorologi¬ 
cal  observations  to  liquid  water  content  profiles.  It  is  designed  for 
computer  applications  and  has  been  applied  to  data  for  the  USSR.  A  deci¬ 
sion  tree  is  used  to  assign  the  liquid  water  content  of  each  atmospheric 
layer.  The  summaries  of  the  analysis  and  procedure  are  presented  in  the 
following  sections. 

2.  Analysis 

The  analysis  consists  of  converting  the  observational  data  for  each 
time  and  place  into  a  liquid  water  content  (LWC)  profile.  The  key  factors 
to  assign  the  liquid  water  content  are  the  temperature  profile,  the  weather 
condition  at  the  surface,  the  present  convective  activity  and  the  precipi¬ 
tation  condition  at  each  level.  In  orde"  to  get  the  LWC  profile,  four 
temperature  regions  (T  <  -30°C,  -30°C  <  T  <_  15°C,  -15°C  <  T  <_  0°C,  T  >  0°C) 
are  used  to  separate  the  precipitation  types,  which  are  combined  with  the 
surface  weather  conditions  (moderate,  heavy  or  light  precipitation)  and 

the  cloud  convective  activity  to  vary  the  liquid  water  content  values.  In 

3 

the  first  region,  ice  water  contents  vary  from  0.1  to  u.2  g/m  depending 
on  the  weather  and  convection.  In  the  second  region,  snow  water  contents 
vary  from  0.4  to  1.2  g/m"*  and  ice  water  contents  are  assigned  from  0.1  to 

3 

0.15  g/m  for  the  layers  without  weather  and  convection.  In  the  third 
region,  snow  water  contents  vary  from  0  to  1.2  g/m'*.  In  the  fourth  region, 
rain  water  contents  vary  from  0  to  0.5  g/m*.  Flag  1  indicates  precipita¬ 
tion  or  convection  without  saturation  in  the  sounding;  the  analysis  then 
has  to  readjust  the  relative  humidity  or  dew  point  depression.  Flag  2 
indicates  moderate  or  heavy  precipitation  without  saturation  at  specific 
layers.  Any  unsaturated  layers  between  saturated  layers  are  also  flagged 
for  analysis  (see  flow  diagram). 
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Procedure 


The  procedure  to  assign  the  LWC  profile  is  as  follows: 

(1)  check  and  correct  observation  data; 

(2)  match  surface  and  radiosonde  observation  for  time,  date 
and  station; 

(3)  categorize  surface  weather; 

a.  precipitation  -  based  on  current  weather  (ww)  in 
surface  observations 

1.  moderate  or  heavy 

ww  (current  weather  code) 

62,  63,  65,  65,  67,  69,  72,  73,  74,  75,  81, 

82,  84,  86,  88,  90,  92,  94,  97,  99 

2.  light 

ww  =  16,  60,  61,  66,  68,  70,  71,  76,  77,  78, 

79,  80,  83,  85,  87,  89,  91,  93,  95,  96,  50,  51, 

52,  S3,  54,  55,  56,  57,  58,  59 

3.  none 

ww  =  other  code  numbers 

b.  convective  activity 

Cl  (low  cloud  type)  =  2,  3,  9 

ww  =  17,  18,  19,  25,  26,  27,  29,  80,  81,  82,  83,  84, 
85,  86,  87,  88,  89,  90,  91,  92,  93,  94,  95,  96,  97, 
98,  99 

(4)  starting  at  the  top  level,  assign  a  liquid  water  content 
according  to  the  flow  diagram  in  figure  1, 


Row  diagram  for  assigning  the  liquid  water  content  to  cloud  and  precipi¬ 
tation  layers.  The  key  to  the  symbols  in  the  diagram  are  provided  in  the 
last  page  of  the  appendix. 


KEY  TO  THE  FLOW  DIAGRAM 


precipitation  in  level  immediately  above  =  wc  t  0 

TEMPERATURE  <  -30°C 

-30°C  <  TEMP  <  - 1 5°C 

-15°C  <  TEMP  <  0oCv 

MODERATE  OR  HEAVY  PRECIPUTATION  AT  GROUND 

LIGHT  PRECIPITATION  AT  GROUND 

CONVECTIVE  ACTIVITY  PRESENT 

SNOW  IN  LEVEL  IMMEDIATELY  ABOVE 

ANY  LEVEL  RH  >_  100 °i 

PRESENT  WEATHER  CODE  =  64,  65,  74,  75,  82,  97 

PRESENT  WEATHER  CODE  =  67,  81,  84,  86,  88,  90,  92,  94 

=  PRECIPITATION  IS  SNOW,  WC  =  a 

=  PRECIPITATION  IS  ICE,  WC  =  a 

■=  PRECIPITATION  IS  RAIN,  WC  -  a 
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